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ABSTRACT 


The  time  variation  of  the  atmospheric  radiation  regime  is 
investigated  for  a  selected  synoptic  situation.  Elsasser's  graphical 
method  shows  the  diurnal  variation  of  the  upward-directed  long-wave  flux 
to  be  in  phase  with  surface  temperature  changes,  but  no  such  change  of 
downward  flux  was  found.  Clouds  reduce  the  upward  flux  above  the  cloud 
tops  and  Increase  the  downward  flux  below  their  bases.  Long-wave 
radiative  cooling  is  greatest  above  cloud  tops  and  at  upper  surfaces  of 
Inversions.  Ground  surface  temperature  is  shown  to  be  very  important  in 
the  long-wave  regime. 

Short-wave  molecular  scattering  and  surface  reflection  are  evaluated 
for  a  model  atmosphere  by  Chandrasekhar's  exact  method  and  are  adapted 
for  the  synoptic  situation.  Water  vapor  absorption,  ozone  absorption, 
and  the  effects  of  atmospheric  dust  are  approximated.  Surface  reflection, 
water  vapor  absorption,  and  absorption  by  particulates  in  a  polluted 
atmosphere  all  play  dominant  roles,  while  the  effects  of  molecular 
scattering  and  ozone  absorption  are  relatively  minor  in  the  over-all 
energy  regime. 
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DISTRIBUTIONS  OF  ATMOSPHERIC  RADIATIVE  HEATING 
AND  COOLING 


I  INTRODUCTION 

The  present  study  is  an  evaluation  of  the  long-wave  and  short-wave 
radiative  energy  transfers  in  a  synoptic  situation  selected  as  being 
typical.  The  motivation  for  the  study  was  two-fold.  First  is  the  need 
for  non-adiabatic  heating  functions  for  use  in  dynamical-numerical 
atmospheric  models.  Secondly,  the  imminent  availability  of  satellite 
radiation  measurements  makes  it  highly  desirable  to  study  the  relation 
between  the  energy  fluxes  emerging  from  the  top  of  the  atmosphere  and 
the  processes  occurring  within  the  atmosphere. 

The  report  covers  a  case  study  of  the  radiative  energy  transfers 
in  a  vertical  column  over  a  given  station  as  a  function  of  time  during 
the  passage  of  a  weather  system.  The  selection  of  a  station  for  the 
study  and  the  general  synoptic  conditions  are  described  in  Sec.  II  of  the 
report.  In  order  to  facilitate  the  computations  and  the  analyses  of  the 
flux  distributions,  the  synoptic  situation  was  idealized;  the  method  and 
assumptions  used  to  accomplish  this  are  discussed  in  Sec.  III. 

The  studies  of  the  long-wave  regime  and  the  short-wave  regime  were 
conducted  separately  and  are  accordingly  reported  separately  in  Secs.  IV 
and  V.  In  the  long-wave  regime,  only  the  radiative  exchanges  due  to  the 
absorption  and  emission  by  water  vapor  are  considered.  The  fluxes  are 
determined  by  means  of  Elsasser's  radiation  chart*^  and  their  distri¬ 
butions  are  presented  in  a  series  of  time  cross  sections.  The  study  of 
the  short-wave  regime  investigates  the  molecular  scattering  and  surface 


*  All  references  are  listed  at  the  end  of  the  report. 


1 


reflection  of  solar  radiant  energy  by  the  use  of  the  exact  theory  of 
radiative  transfer  as  first  derived  by  Chandrasekhar.®  Water  vapor 
absorption  of  short-wave  radiation  is  examined  by  an  empirical  method 
developed  by  Korb,  £t  al.,*^  and  the  effects  of  haze  and  dust  are 
determined  by  means  of  a  three-layered  atmospheric  model.  The  rate  of 
heating  by  ultraviolet  absorption  of  ozone  is  calculated  on  the  basis 
of  the  latest  information  on  ozone  distribution. 

Concurrent  investigations  under  this  project  are  concerned  with 
the  dynamical  significance  of  non-adlabatlc  heating  and  cooling  distri¬ 
butions  on  the  evolution  of  atmospheric  circulations.  Dynamical- 
numerical  Integration  with  prescribed  fields  of  heating  rates  are  carried 
out  for  this  purpose.  This  w6rk  will  appear  in  a  later  report.  These 
parallel  investigations  should  lead  to  dependent  heating  formulations 
which  will  be  useful  in  further  realization  of  the  Inherent  predicta¬ 
bility  of  atmospheric  motions  and  the  weather. 


2 


II  THE  SYNOPTIC  SIITJATION 


The  selection  of  a  station  for  the  case  study  required  that  the 
station  have  (l)  sufficient  meteorological  data  to  facilitate  the  compu¬ 
tations  of  radiation  fluxes,  (2)  a  cyclone  or  frontal  passage  to  bring 
about  changes  In  the  state  of  the  atmosphere,  (s)  clear-sky  conditions 
before  and  after  the  passage  of  the  disturbance  as  "references,**  and 
(4)  a  location  preferably  on  flat,  continental  terrain. 

The  station  selected  as  meeting  these  r''qulrements  was  Omaha, 
Nebraska,  which  experienced  a  cyclone-frontal  passage  during  the  period 
12-16  March  1957.  At  the  beginning  of  the  period  (12  March  1957)  almost 
the  entire  Midwest,  Including  Omaha,  was  under  a  weak  high-pressure 
system  and,  except  for  a  few  scattered  cirrus,  clear  skies  prevailed 
throughout  most  of  the  area.  On  the  13th  of  March,  a  weak  low  situated 
on  the  lee-slde  of  the  Rocky  Mountains  Initiated  a  regeneration  with  an 
accompanying  frontogenesls .  This  continued  to  develop  on  the  14th,  the 
center  of  the  low  and  the  polar  front  with  Its  associated  weather  passing 
through  Omaha  at  about  0300Z  on  this  date.  This  was  followed  on  the  15th 
and  16th  by  the  outbreak  of  cold  continental  polar  air  behind  the  front 
and  a  bulld-up  of  a  high-pressure  cell  with  clear  skies  over  Omaha  and 
much  of  the  Midwest.  'Figure  1  shows  the  synoptic  situation  over  the 
United  States  at  1830Z,  14  March  1957. 

The  time-section  analysis  of  the  synoptic  situation,  as  well  as  the 
flux  computations  and  the  analyses  of  the  radiation  fields,  was  based  on 
the  usual  meteorological  data  (pressure,  temperature,  humidity)  as  ‘ 
published  In  the  U.S.  Weather  Bureau's  Northern  Hemisphere  Data  Tabulation 
and  on  synoptic  data  (WBAN  lOA  and  b)  obtained  from  the  National  Weather 
Records  Center.  Lack  of  data  on  the  moisture  distribution  above  400  mb 
made  It  necessary  to  extrapolate  the  moisture  curve  to  the  tropopause 
level  In  all  cases.  The  procedure  used  here  was  to  extend  the  dewpoint 
sounding  upward  from  the  point  of  Its  termination  with  a  curve  that 
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SYNOPTIC  SITUATION  OVER  U.S.,  1830Z,  14  MARCH  1957 
From  USWB  Daily  Weather  Map  Series,  U.S.  Dept,  of  Commerce 


paralleled  the  mean  moisture  distribution  at  the  latitude  for  Omaha  as 
given  by  London.* 

Insufficient  data  on  the  heights  of  the  cloud  bases  and  tops  also 
posed  a  serious  obstacle  in  the  determination  of  the  radiation  fluxes. 
Some  information  on  cloud-base  heights  was  available  from  the  synoptic 
reports,  but  values  of  cloud- top  heights  were  completely  lacking. 
Estimates  were  made  by  using  the  synoptic  reports  and  the  cloud  analysis 
criteria  employed  by  Wilk*  as  guides  and  subjectively  analyzing  the 
meteorological  soundings. 
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Ill  MODELING  OF  THE  SYNOPTIC  SITUATION 


The  usual  synoptic  time  cross-sectional  analysis  of  the  meteorological 
parameters  (i.e.,  pressure,  temperature,  moisture  plus  clouds,  frontal 
surface  and  weather)  presented  much  too  complex  a  picture  of  the  changing 
atmospheric  state  for  computing  radiative  fluxes.  In  order  to  make  it 
more  tractable  for  our  study,  the  original  time  cross-sectional  analysis 
was  simplified  by  reproducing  the  synoptic  situation  in  a  schematic 
representation  as  shown  in  Fig.  2.  In  the  model,  the  vertical  distributions 
of  pressure,  temperature,  and  moisture  as  reported  in  the  observational 
data  were  faithfully  retained,  for  the  most  part,  but  the  field  of 
each  parameter  was  smoothed  of  some  detail.  An  observation  of  blowing 
dust  reported  during  one  period  was  omitted  because  it  added  complications 
requiring  study  that  was  beyond  the  scope  of  the  present  effort. 

The  most  significant  modification  of  the  synoptic  conditions 
was  made  with  respect  to  the  cloud  distributions.  Besides  the  previously 
mentioned  lack  of  data  on  the  heights  of  cloud  tops,  the  complexity 
introduced  into  the  radiation  field  by  the  physical  parameters  of  the 
clouds  compelled  gross  simplifications  to  be  made  as  follows:  (l)  all 
clouds  were  assumed  to  be  black  bodies  and  only  radiation  from  their 
tops  and  bases  was  considered,  (2)  all  amounts  of  cloudiness  were 
considered  as  moving  over  the  station  as  one  smooth-surfaced  cloud 
body,  (3)  the  period  for  the  radiation  field  to  go  from  the  clear-sky 
condition  to  cloudy  condition  was  assumed  to  be  two  hours  for  all  amounts 
of  cloudiness,  (4)  a  discontinuity  in  the  flux  distribution  was  considered 
to  occur  right  at  the  leading  and  trailing  edges  of  the  cloud  layers, 

(5)  heavy  continuous  precipitation  was  also  assumed  to  be  black  body  and 
was  considered  together  with  the  accompanying  cloud  as  such,  and  (6) 
multi-layered  clouds  were  assumed  continuous  in  the  vertical  such  as  to 
form  a  single  layer.  In  the  case  of  the  partial  cloudiness  the  fluxes 
were  determined  by  weighting  the  fluxes  according  to  the  following  formula: 
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FIG.  2 

SCHEMATIC  REPRESENTATION  OF  THE  SYNOPTIC  SITUATION  AT 
Thin  solid  lines  are  isotherms  in  "C,  dash-dot  iines  are  i! 
in  g/kg,  front  and  inversions  are  represented  by 
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FIG.  2 

REPRESENTATION  OF  THE  SYNOPTIC  SITUATION  AT  OMAHA  DURING  12-16  MARCH  1957 
id  lines  are  Isotherms  in  “C,  dash-dot  lines  are  isopleths  of  constant  mixing  ratio 
in  g/kg,  front  and  inversions  are  represented  by  heavy  solid  lines. 


Fp  =  F^(l  -  n)  +  Fq(n) 

where  Pp  is  the  weighted  flux  value  for  partial  cloud  cover  condition, 

Ff,  is  the  flux  value  for  clear  sky  condition,  F  ,  is  the  flux  value  for 
overcast  conditions,  N  is  the  amount  of  cloud  cover  in  tenths  of  the  sky. 
N  =  1.0  indicates  overcast  conditions,  and  N  =  0  indicates  no  cloud 
conditions. 

There  is  no  measurement  on  a  regular  routine  basis  of  the  soil 
surface  temperature  which  is  needed  for  determining  the  flux  of  radiant 
energy  emitted  by  the  earth's  surface.  Thus  in  the  time  section,  the 
ground  surface  temperature  was  assumed  to  be  equivalent  to  the  reported 
surface  air  temperature  of  the  meteorological  observations .  The  effects 
and  consequences  of  this  assumption  are  discussed  in  a  later  section. 
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IV  THE  LONG-WAVE  TERRESTRIAL  RADIATION  REGIME 


A.  GENERAL 

The  primary  interests  in  the  study  of  the  long-wave  terrestrial 
radiation  were  the  determination  of  the  long-wave  flux  distribution  and  its 
variations  within  the  troposphere  during  changing  synoptic  conditions, 
the  examination  of  the  physical  processes  responsible  for  bringing 
about  the  variations,  and  the  determination  of  the  heating  or  cooling 
occurring  in  the  troposphere  resulting  from  long-wave  radiative  exchanges. 
Another  aspect  of  the  study  was  to  relate  the  emergent  radiant  fluxes 
from  the  top  of  the  atmosphere,  as  might  be  detected  by  a  satellite,  with 
the  radiation  processes  within  the  atmosphere. 

The  atmospheric  constituents  that  are  Important  for  the  absorption 
of  radiant  energy  within  the  spectral  range  of  terrestrial  radiation  are 
water  vapor,  carbon  dioxide,  and  ozone.  Although  at  any  level  (especially 
the  higher  levels)  of  the  troposphere,  the  carbon  dioxide  flux  may  be  a 
significant  percentage  of  either  the  total  upward  or  downward  long-wave 
flux;  its  net  contribution  to  long-wave  radiative  temperature  changes  in 
the  troposphere  is  generally  very  minor  compared  to  that  of  water  vapor. 

For  our  purposes,  it  was  considered  insignificant.  Similarly  the  effects 
of  tropospheric  ozone  were  not  Included  in  the  present  study  since  the 
small  concentration  of  this  gas  makes  its  role  in  the  tropospheric  radiative 
cooling  or  heating  processes  inappreciable.  Thus,  the  present  study  of 
terrestrial  radiation  considered  only  the  radiative  exchanges  due  to  the 
absorption  and  emission  by  water  vapor  which,  in  addition  to  its  important 
absorptive  properties,  is  distributed  within  the  troposphere  in  greater 
amounts  than  the  other  two  absorbing  gases. 

B.  METHOD  OF  COMPUTATION 

Elsasser's®  graphical  method  was  employed  in  the  computations  of  the 
upward  flux,  F^,  and  the  downward  flux,  F^j,  and  from  these  two  the  net 
flux,  Fj^,  and  the  flux  divergence  were  determined.  While  it  is  understood 
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that  later,  improved  methods  of  computing  fluxes  of  long-wave’ 
radiation  are  available,  the  use  of  Elsasser's  earlier  technique  is 
justified  on  the  basis  of  (l)  convenience,  and  (2)  the  fact  that  the 
values  obtained  are  thought  to  be  sufficiently  accurate  to  give  a  valid 
pattern  and  the  essential  changes  of  the  distribution  of  the  radiation 
parameters . 

The  optical  path  length,  which  was  used  in  computing  the  fluxes  from 
the  radiation  charts,  was  corrected  for  the  pressure-dependence  of  line¬ 
broadening  in  the  water-vapor  spectrum  according  to  London^ .  In  its 
final  form,  the  corrected  optical  path  length,  u,  is  defined  as: 


where  ^  is  the  acceleration  due  to  gravity,  w  is  the  mixing  ratio, 

Pi  and  Pa  are  the  pressures  at  lower  and  upper  boundaries,  respectively, 

of  the  atmospheric  layer,  and  p  is  the  standard  pressure,  1000  mb. 

s 

Flux  computations  were  made  to  the  tropopause  (250-to-300-mb  region) 
using  the  mandatory  levels  of  the  radiosonde  reports  as  reference  levels. 
Computations  were  made  at  additional  levels  at  those  times  when  an 
inversion,  a  frontal  surface,  or  a  cloud  was  present.  Since  the 
meteorological  data  were  available  only  at  12-hour  intervals,  supplementary 
soundings  were  composed  from  the  synoptic  time  section  for  6-hour  periods. 
Flux  computations  were  then  made  for  the  shorter  periods  in  order  to 
aid  in  the  analyses  of  the  radiation  flux  fields. 

C .  UPWARD  FLUX 

Figure  3  represents  the  time-sectional  analysis  of  the  upward  flux, 

F^,  of  terrestrial  radiation.  It  displays  definite  patterns  of 
distribution  which  show  that  F^  generally  decreases  with  height  and  that 
the  largest  decreases  take  place  in  the  lower  troposphere  where  the 
water-vapor  concentration  is  the  greatest. 


10 


THE  UPWARD  LONG  WAVE  RADIATIVE  FLUX 
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FIG.  3 

DISTRIBUTION  OF  THE  UPWARD  LONG-WAVE  RADIATIVE  FLUX  . 
Units:  ly/3  hr.  The  dashed  lines  represent  flux  values  ini 

given  by  the  solid  lines. 
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The  most  prominent  features  of  the  time  section  are  the  marked 
diurnal  variation  of  and  the  effect  of  clouds  on  the  vertical  distri¬ 
bution  of  Fy.  The  diurnal  variation  of  F^  is  in  phase  with  the  ground 
surface  temperature  such  that  the  maximum  at  any  one  level  occurs  at  the 
time  of  the  midday  maximum  of  the  ground  surface  temperature  and  the 
minimum  occurs  at  the  time  of  the  early-morning  minimum  of  the  surface 
temperature.  This  "in  phase"  time  variation  of  the  surface  temperature 
an^.  the  vertical  distribution  of  the  F^  is  the  natural  consequence  of 
the  earth's  surface  radiating  as  a  black  body.  It  is  evident  throughout 
the  cloudless  troposphere,  but  the  absorption  of  terrestrial  radiation 
by  water  vapor  somewhat  dampens  the  relationship  at  the  upper  levels. 

The  effect  of  clouds  is  to  decrease  F^,  from  clear  sky  conditions, 
in  the  region  of  the  atmosphere  above  the  cloud  tops.  The  largest 
gradients  in  the  distribution  of  F^  occur  above  the  edges  of  the  cloud 
layers  and  the  time  rate  of  change  is  greatest  at  those  levels  nearest 
to  the  cloud  tops.  The  uniform  distribution  of  F^  above  the  clouds 
reflects  the  uniformity  of  the  cloud- top  temperatures  which  were  assumed 
to  be  the  same  as  the  free-air  temperatures.  The  black-body  assumption 
for  the  cloud  surfaces  results  in  the  vertical  distribution  of  Fy  ap¬ 
proaching  clear-sky  conditions  with  decreasing  height  of  the  cloud  tops. 

With  respect  to  the  cloudless  atmosphere,  there  is  some  change  in 
the  magnitudes  and  in  the  distribution  of  F^  between  the  warm,  moist  air 
ahead  of  the  front  and  the  cold,  dry  polar  air  behind  it.  In  the  former, 
the  flux  values  are  generally  larger  throughout  the  air  mass  because  of 
the  higher  radiating  temperatures  [60  to  90  ly  (3  hr)“*  as  compared  to 
50  to  70  ly  (3  hr)“*  behind  the  front],  and  the  decrease  of  F^  with  height 
is  greater  due  to  increased  absorption  brought  about  by  the  larger  water- 
vapor  content.  The  polar  air  mass  following  the  front  exhibits  the  same 
diurnal  variation  pattern  of  F^  as  in  the  pre-frontal  air  but  its  drier 
nature  results  in  greater  transparency  to  terrestrial  radiation.  The 
decrease  of  F^  through  a  column  from  the  earth's  surface  to  the  tropopause 
is  of  the  order  of  7  to  15  ly  (3  hr)"^  in  the  polar  air  and 
15  to  25  ly  (3  hr)”*  in  the  pre-frontal  air. 


12 


There  is  no  consistent  trend  in  the  distribution  through  inversions. 
Results  of  the  analysis  show  both  increases  and  decreases  averaging  about 
2.0  ly  (3  hr)"*  from  the  bases  to  the  tops  of  inversions,  depending  on  the 
strength  (temperature  gradient)  and  the  moisture  distribution  through  the 
stable  layer. 

D.  UPWARD  FLUX  AT  THE  TROPOPAUSE 

The  time  variation  of  Fu  passing  through  the  tropopause  level  is  shown 
by  the  curve  in  the  upper  section  of  Fig.  3.  The  curve  essentially  repre¬ 
sents  the  time  variation  of  F^  emerging  from  the  top  of  the  atmosphere, 
since  the  inclusion  of  radiative  exchanges  by  the  small  amount  of  water 
vapor  and  by  the  other  absorbing  gases  would  not  significantly  alter  the 
configuration  of  the  curve,  although  the  magnitudes  may  be  changed.  For 
instance.  Hales  and  Zdunkowski^  have  computed  the  radiation  fluxes  of 
water  vapor,  carbon  dioxide,  and  ozone  leaving  the  earth's  atmosphere  for 
different  types  of  air  masses.  While  no  direct  comparison  could  be  made 
because  their  computations  were  made  for  different  atmospheric  conditions, 
their  values  of  radiative  flux  at  the  upper  limit  of  the  atmosphere  are 
less,  by  factors  of  1.5  to  2.0,  than  the  values  of  the  upward  flux  at  the 
tropopause  obtained  in  the  present  study.  On  the  other  hand,  their  results 
show  that  the  contribution  to  the  total  flux  leaving  the  top  of  the  atmos¬ 
phere  by  water  vapor,  carbon  dioxide,  and  ozone  are  in  about  the  same  ratio 
for  all  air  masses  considered. 

The  relationship  between  the  upward  flux  at  the  tropopause  and  the 
radiative  processes  within  the  atmosphere  may  be  seen  from  the  upper  curve 
and  the  time  cross  section  of  Fig.  3.  The  largest  values  of  the  emergent 
flux  are  found  in  the  cloudless  atmosphere  at  the  time  of  the  midday 
maximum  of  the  ground  surface  temperature.  The  diurnal  variation  of  the 
flux  still  oscillates  in  phase  with  the  varistion  of  the  ground  surface 
temperature  but  with  less  magnitude  of  change  than  in  the  lower  tropo¬ 
spheric  levels.  The  outstanding  feature  is  the  sharp  decrease  of  the 
outward  flux  from  clear-sky  conditions  at  the  periphery  of  the  cloud 
layers,  with  a  more  or  less  uniform  distribution  in  the  area  above  the 
cloud  top  proper.  The  decrease  becomes  greater  with  increasing  height 
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of  the  cloud  tops  because  of  the  lower  cloud-top  temperatures  found  at 
the  higher  levels. 

E .  DOWNWARD  PLUX 

The  changes  In  the  distribution  of  the  downward  flux,  F^j,  are  not 
as  pronounced  as  the  ones  depicted  by  the  distribution  patterns  of  Fy. 
However,  some  general  remarks  can  be  made  from  the  time-sectional  analysis 
of  the  Fjj  distribution  represented  in  Fig.  4. 

An  examination  of  the  field  shows,  first  of  all,  that  F^  increases 
from  the  tropopause  to  the  earth's  surface.  There  is  no  regular  diurnal 
variation  of  F^  and  the  variations  in  the  distribution  of  F^j  throughout 
the  day  are  small  compared  to  the  diurnal  changes  of  F^.  Instead  of  the 
temperature  dependence  characterized  by  the  distribution  of  F^,  the  Iso- 
pleths  of  F^  parallel  the  mixing  ratio  Isopleths  in  the  region  below 
5  km,  indicating  the  close  relationship  of  F^j  and  the  moisture  distribution. 

At  the  very  high  levels,  the  isopleths  of  F^  appear  to  parallel  the 
tropopause.  This  reflects,  at  least  in  part,  the  restrictions  of  the 
computational  method,  namely,  not  taking  into  account  the  water  vapor 
above  the  tropopause.  On  the  other  hand,  the  small  amount  of  water  vapor 

in  the  lower  stratosphere  would  not  change  the  flux  values  by  an  ap- • 

preclable  amount.  Hales  and  Zdunkowskl^  found  that  changing  the  temper¬ 
ature  and  moisture  distribution  at  high  altitudes  in  several  different 
ways  did  not  result  in  any  significant  changes  of  the  water-vapor  flux 
density. 

The  presence  of  a  cloud  produces  a  marked  change  in  the  vertical 
distribution  of  F^j.  Its  effect  is  to  increase,  F^j,  relative  to  clear  sky 
condition,  below  the  cloud  layer,  with  the  largest  increase  occurring  just 
below  the  cloud  base.  With  increasing  optical  path  length  from  the  cloud 

base,  the  cloud  effect  on  the  vertical  distribution  of  F^  is  reduced  and, 

in  the  cases  of  partial  cloudiness  at  1500Z,  12  March  and  0300Z,  13  March, 
in  the  time  section,  it  is  completely  damped  out  at  about  5  km  below  the 
cloud  layers.  The  flux  that  reaches  the  earth's  surface  is  affected  more 
by  clouds  with  low  bases  than  by  those  with  bases  at  higher  elevations. 
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THE  DOWNWARD  LONG  WAVE  RADIATIVE  FLUX 
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)N  OF  THE  DOWNWARD  LONG-WAVE  RADIATIVE  FLUX  AT  OMAHA  DURING  12' 

Units:  ly/3  hr 


16  MARCH  1957 


Analogous  to  the  change  In  the  upward  flux  distribution  above  cloud  layers, 
the  largest  change  in  the  distribution  of  Fj  occurs  at  the  leading  or 
trailing  edges  of  the  cloud  bases. 

The  larger  .Mounts  of  water  vapor  in  the  lower  troposphere  (<  5  kin) 
result  in  the  rate  of  change  of  in  the  vertical  direction  being  greater 
there  than  in  the  upper  troposphere.  This  effect  of  Increased  water  vapor 
is  also  shown  by  the  difference  in  the  gradients  between  the  air  nasses. 

In  the  nolst  air  ahead  of  the  front,  the  flux  increases  by  approzloately 
30  ly  (3  hr)~'  from  5  km  to  the  ground,  while  in  the  drier'  air  behind 
it,  the  Increase  between  the  same  two  levels  is  a  little  over  20  ly 
(3  hr)”*.  The  difference  between  the  vertical  gradients  of  Fj  of  the  two 
air  Basses  is  not  as  great  for  the  tropospheric  region  from  the  tropopause 
to  the  b-kn  level. 

The  computed  values  of  F^j  in  the  clear  sky  show  decreases  of  0.5  to 
3.0  ly  (3  hr)‘“  from  the  bases  to  the  tops  of  inversions.  The  snail 
changes  of  F()  with  tine  mask  the  decrease  of  the  flux  through  stable 
layers,  and  thus  the  decrease  does  not  appear  in  the  time  section. 

F.  DOWNWARD  FLUX  AT  THE  GROUND 

The  curve  shown  in  the  upper  section  of  Fig.  4  represents  the  time 
variation  of  Fa  reaching  the  ground.  In  the  cloudless  sky,  the  change  is 
not  too  large  over  relatively  long  periods;  for  Instance,  between  1500Z, 

12  March  and  1500Z,  13  March,  the  flux  values  for  the  period  vary 
5:5  ly  (3  hr)"*  from  the  value  at  the  start  of  the  period.  There  is 
partial  cloudiness  at  1500Z,  12  March  and  0300Z,  13  March,  but,  as 
mentioned  before,  the  Increase  of  F^  from  the  cloud  bases  is  dampened  by 
the  presence  of  water  vapor  in  the  lower  atmosphere  and  the  effect  of 
clouds  on  F^s  is  not  apparent  at  the  ground.  The  most  significant  change 
is  the  large  increase  of  Fj  brought  about  by  the  cloud  layer  during  the 
period  ISOOZ,  13  March  to  0900Z,  14  March.  The  gradual  change  of  F,j 
reaching  the  earth's  surface  at  the  beginning  of  the  period  is  due  to 
(l)  the  gradual  lowering  of  the  cloud  base,  thus  a  gradual  change  of  the 
radiating  temperature  of  the  cloud  base,  and  (2)  the  decrease  of  optical 
path  length  between  the  cloud  base  and  the  earth's  surface  over  a 
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relatively  extended  period.  A  more  abrupt  change  from  overcast  to  clear 
sky,  and  then  to  overcast  conditions,  Is  displayed  by  the  sharp  decrease 
and  Increase  of  the  curve  during  0300Z  to  0900Z,  14  March. 

G.  FLUX  DIVERGENCE 

The  rate  of  cooling  was  determined  according  to  the  following 
expression  (Elsasser®): 


where  g  Is  the  acceleration  due  to  gravity,  Cp  Is  the  specific  heat  of 
air  at  constant  pressure,  and  Pj  represent  the  net  flux  and  the 
pressure,  respectively,  at  the  lower  boundary  of  the  atmospheric  layer, 
and  Fjj^  and  Pj  are  the  net  flux  and  pressure,  respectively,  at  the  upper 
boundary  of  the  atmospheric  layer.  The  analysis  of  the  rate  of  cooling 
for  the  synoptic  situation  under  Investigation  Is  shown  In  Fig.  5. 

The  field  of  flux  divergence  or  rate  of  cooling  displays  more 
cellular  patterns  than  does  the  upward  or  downward  flux  field.  It  Is 
difficult  to  establish  any  regular  variations  In  the  cooling  patterns 
or  In  the  vertical  distribution  of  radiative  cooling  In  the  clear  sky. 

In  a  very  general  sense,  there  appears  to  be  relatively  large  cooling  In 
the  cloudless  atmosphere  below  2.5  km;  above  this  level  the  rate  of 
cooling  decreases  to  about  5  km  and  then  Increases  to  a  point  below  the 
tropopause  level . 

The  largest  rates  of  cooling  are  found  above  the  tops  of  cloud  layers 
and  at  or  above  the  upper  boundaries  of  Inversions  and  fronts.  Even  with 
partial  cloudiness,  the  cooling  above  cloud  layers  Is  Increased  as  compared 
with  either  cloudless  conditions  or  with  the  cooling  rates  just  below  the 
cloud  bases.  The  computed  values  of  the  radiative  cooling  rates  above 
cloud  tops  range  from  0.08®C  (3  hr)”*  for  0.6  cloud  cover  at  1500Z, 

12  March,  to  a  maximum  of  0.40*C  (3  hr)"*  for  the  overcast  at  1500Z, 

14  March. 
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THE  NET  LOSS  OF  ENERGY  BY  THE  TROPOSHERE 
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FIG.  5 

DISTRIBUTION  OF  LONG-WAVF  FLUX  DIVERGENCE  (RATE  OF  CC 
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FIG.  5 

F  LONG-WAVE  FLUX  DIVERGENCE  (RATE  OF  COOLING)  AT  OMAHA  DURING  12-16  MARCH  1957 

Units;  *C/3  hr 


Below  the  cloud  bases,  less  cooling  takes  place  through  fairly  large 
depths  (2  to  4  km)  of  the  troposphere  below  the  clouds.  The  overcast 
cover  and  the  low  level  of  the  cloud  base  (higher  radiating  temperature) 
result  in  a  region  of  radiative  heating  below  the  cloud  layer  during  the 
period  1500Z,  13  March  to  0300Z,  14  March.  A  small  region  of  slight 
heating  is  also  Indicated  below  the  low  cloud  base  at  0900Z,  14  March. 
Below  the  clouds,  the  computed  temperature  changes  range  from  a  cooling 
of  O.OS^C  (3  hr)"*  to  a  heating  of  0.14®C  (3  hr)"*. 

The  centers  of  extreme  cooling  above  and  near  the  leading  edges  of 
clouds  and  the  sharp  gradient  and  discontinuity  at  the  leading  and 
trailing  edge  of  the  cloud  base  ((jl 500 Z;,  13  March;  0300Z,  14  March)  are 
due  to  the  assumption  of  a  discontinuity  in  the  flux  distribution  at  the 
periphery  of  the  cloud  layers.  However,  the  general  pattern  of  cooling 
above  the  cloud  layers  and  less  cooling  or  heating  below  the  cloud  bases 
are  in  agreement  with  the  results  of  Alder.®  He  found  maximum  radiative 
cooling  to  occur  2000  to  3000  ft  above  cloud  tops,  and  maximum  heating 
below  cloud  layers  to  occur  at  the  bases. 

The  largest  rate  of  cooling  [0.50*0  (3  hr)”*]  is  found  at  the  upper 
boundary  of  the  frontal  surface.  This  is  due  largely  to  the  large  amount 
of  water  vapor  in  the  frontal  zone,  with  a  rapid  decrease  of  water  vapor 
above.  ^  The  rate  of  cooling  decreases  quite  sharply  through  the  frontal 
layer  and  near  the  ground.  The  lower  boundary  exhibits  a  slight  heating. 
In  a  like  manner,  increased  cooling  is  seen  to  take  place  at  or  above  the 
tops  of  Inversion  layers,  with  less  cooling  at  the  bases.  The  values 
range  from  0.20  to  0.36“C  (3  hr)”*  at  the  top,  to  0.06  to  0.24®C  (3  hr)”* 
at  the  bases  of  the  inversions. 

Computations  were  made  to  determine  the  rate  of  cooling  within 
clouds  by  employing  Eq.  (l),  in  which  a  uniform  rate  of  cooling  is  assumed 
to  take  place  throughout  the  whole  depth  of  the  cloud  layer.  No  study  of 
the  cooling  within  a  cloud  was  made  other  than  to  get  an  idea  of  the  order 
of  magnitude.  Hoffer®,  using  the  same  expression,  found  that  the  rate  of 
cooling  within  a  cloud  was  approximately  a  hyperbolic  function  of  the 
pressure  thickness  of  the  cloud  layer  and  that  if  two  clouds  were  of  the 
same  pressure  thickness,  then  the  cloud  with  the  base  at  a  higher  pressure 
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cooled  more  rapidly.  The  computed  values  of  cooling  are  shown  within 
each  cloud  layer  in  the  time  section  of  Fig.  5. 

H.  LOSS  OF  ENERGY  BY  THE  TROPOSPHERE 

The  curve  above  the  time  section  of  Fig.  5  shows  the  time  variation 
of  the  net  energy  loss  by  long-wave  radiative  processes  from  an  atmos¬ 
pheric  column  extending  from  the  earth's  surface  to  the  tropopause.  The 
positive  values  of  the  curve  show  that  although  there  may  be  sub-synoptic 
regions  of  long-wave  radiative  cooling  and  heating,  the  net  effect  of  the 
long-wave  radiative  processes  is  the  loss  of  energy  by  the  troposphere. 

The  loss  is  relatively  constant  throughout  the  day  in  the  cloudless 
atmosphere.  The  large  variations  of  the  energy  loss  are  introduced  by 
the  presence  of  overcast  cloud  cover. 

The  minimum  amount  of  energy  loss  occurs  during  the  period  from 
1500Z,  13  March  to  about  OlOOZ,  14  March.  During  the  period,  the  net 
flux,  Fjj,  at  the  tropopause  remains  more  or  less  constant,  so  the  vari¬ 
ations  of  energy  loss  are  dependent  on  what  takes  place  at  and  below  the 
cloud  layers.  Between  1500Z  and  2100Z  on  the  13th,  the  cloud  base  lowers, 
thus  increasing  its  radiating  temperature  and  F^j,  which  reaches  the  ground. 
At  the  same  time,  the  ground  surface  temperature  increases  with  an  ac¬ 
companying  Increase  of  Fu(see  Fig.  2)  Since  Fu  at  the  ground  Increases 
more  than  Fjj,  F^  is  changed  so  that  there  is  a  net  gain  of  energy  by  the 
atmosphere  at  the  earth-atmosphere  It"  i  Cace,  and  a  decrease  of  the  net 
energy  loss  from  the  troposphere.  F^om  2100Z,  13  March  to  about  OlOOZ 
14  March,  the  cloud  base  remains  at  a  constant  level,  as  does  the  cloud 
top,  but  the  decrease  of  the  ground  s  -face  temperature  results  in  Fu 
and  Fjj  at  the  earth's  surface  being  decreased.  Therefore,  during  this 
period,  there  is  less  net  energy  transferred  to  the  atmosphere  from  the 
ground,  or  an  increase  in  the  net  loss  of  energy  by  the  troposphere. 

The  maximum  of  energy  loss  by  the  troposphere  occurs  from  0900  to 
1300Z,  14  March.  At  0900Z,  the  low  height  of  the  cloud  top  results  in 
a  relatively  high  radiating  temperature  and  a  fairly  large  value  of  F^ 
and  Fj^  at  the  tropopause.  Below  the  cloud  layer  the  temperature  of  the 
cloud  base  is  higher  than  that  of  the  ground  surface  and  results  in  Fjj 
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at  the  earth's  surface  being  directed  from  the  atmosphere  to  the  ground 
instead  of  vice-versa,  as  is  generally  the  case.  This  loss  of  energy  at 
0900Z  is  comparatively  large.  The  black-body  assumption  for  the  cloud- 
precipitation  combination  results  in  the  lowering  of  the  black-body 
surface  (cloud  base)  to  the  ground  during  the  per^''d  from  about  lOOOZ  to 
1300Z.  The  consequent  change  in  the  net  energy  loss  by  the  troposphere 
is  a  slight  decrease  from  the  amount  at  0900Z. 

I .  EFFECTS  OF  GROUND  SURFACE  TEMPERATURE 

It  was  assumed  in  the  present  study  that  the  radiating  temperature 
of  thcs  ground  surface  was  equivalent  to  the  reported  surface  air  temper¬ 
ature  given  in  the  regular  meteorological  observation.  This  has  two 
consequences.  First,  it  is  known  from  observations  (Geiger 7> Lettau  and 
Davidson?) that  the  temperature  of  the  earth's  surface  can  differ  appreci¬ 
ably  from  the  temperature  of  the  air  layer  next  to  it.  Since  the  flux  of 
radiant  energy  emitted  by  the  earth's  surface  is  proportional  to  the 
fourth  power  of  its  absolute  temperature,  the  amount  of  energy  computed 
by  the  application  of  the  assumption  above  may  result  in  a  significant 
error.  Computations  of  fluxes  were  made  to  determine  the  probable  error 
that  may  arise  from  neglecting  any  air-ground  temperature  differences. 

The  changes  in  the  vertical  distributions  of  F^  and  F^^  for  5“C  differences 
of  the  ground  temperature  from  that  assumed  in  the  study  are  shown  in 
Fig.  6.  The  change  in  F^  and  F^^  due  to  different  radiating  temperatures 
of  the  ground  is  greatest  at  the  earth-atmosphere  interface,  decreases 
slightly  with  height  in  the  lowest  air  iayers7  and  is  more  or  less  con¬ 
stant  throughout  the  tropospheric  levels  above.  The  curves  show  that  if 
the  ground  temperature  is  greater  than  the  assumed  temperature,  the 
radiant  flux  emitted  by  the  earth's  surface  is  increased  by  approximately 
eight  percent  for  each  S^C  and  is  decreased  by  nearly  the  same  amount  when 
the  ground  temperature  is  less  than  the  assumed  temperature.  These  vari¬ 
ations  of  Fy  result  in  F^^  at  the  interface  being  increased  or  decreased 
by  about  twenty-five  percent  for  each  5“C  increment  of  temperature 
difference . 
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FIG.  6 

THE  EFFECT  OF  THE  GROUND  SURFACE  TEMPERATURE  ON  THE  UPWARD  AND  NET  LONG-WAVE  FLUX  DISTRIBUTIONS 


At  the  middle  and  upper  tropospheric  levels,  Including  the  tropopause, 
and  are  changed  In  the  same  sense  as  mentipned  above,  but  the  changes 
are  only  4^6  to  5^  for  each  5®C  temperature  difference. 

The  second  consequence  introduced  by  assuming  the  ground  surface 
temperature  to  be  the  same  as  the  reported  surface  air  temperature  is  that 
of  a  phase  displacement  of  the  diurnal  oscillation.  Generally,  the  phase 
difference  between  the  diurnal  variations  of  the  two  temperatures  is  of 
the  order  of  only  an  hour  or  so.  This  corresponds  to  a  shift  of  the 
patterns  by  a  like  amount  in  the  time-section  representation,  while  the 
distribution  patterns  themselves  are  otherwise  unaffected  by  the  shift. 

J .  SUMMARY 

The  case  study  of  the  long-wave  terrestrial  radiation  regime  shows 
that  the  distribution  of  the  upward  long-wave  flux  in  the  clear  atmosphere 
decreases  with  height  and  has  a  time  variation  that  is  in  phase  with  the 
diurnal  oscillation  of  the  ground  surface  temperature.  The  downward  flux 
also,  during  cloudless  conditions,  increases  from  the  tropopause  to  the 
earth's  surface  but  does  not  exhibit  any  regular  diurnal  variations  and, 
in  addition,  its  changes  with  time  are  small  compared  to  those  of  the 
upward  flux.  The  most  important  feature  for  producing  pronounced  changes 
in  the  distributions  of  the  fluxes  is  the  presence  of  clouds.  The  effect 
of  clouds  is  to  reduce  the  upward  flux,  from  clear  sky  conditions,  in  the 
region  above  the  cloud  tops.  The  black-body  assumption  for  the  cloud 
surfaces  results  in  the  distribution  of  the  upward  flux  above  the  cloud 
tops  to  approach  that  of  cloudless  conditions  with  decreasing  height  of 
the  cloud  tops.  The  downward  flux  is  increased  by  clouds  in  the  region 
below  the  cloud  bases.  The  flux  reaching  the  ground  is  decreased  with 
increasing  cloud  height,  because  of  lower  cloud-base  temperature  and  more 
absorption  by  water  vapor.  The  time-varying  field  of  flux  divergence  is 
of  a  more  cellular  nature  than  either  the  upward  or  downward- flux  field. 

In  the  cloudless  atmosphere,  the  rate  of  cooling  is  relatively  large  in 
the  air  layers  below  about  2.5  km,  decreasing  above  that  level  to  the 
middle  troposphere,  and  then  increasing  to  the  tropopause.  Maxima  of 
cooling  occur  above  cloud  layers  and  at  or  above  the  upper  boundaries  of 
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inversions  and  fronts.  Less  cooling  or  heating,  depending  on  the  synoptic 
conditions,  takes  place  in  the  region  below  cloud  layers  and  at  the  lower 
boundaries  of  inversions  and  fronts.  The  contribution  of  the  long-wave 
radiative  exchanges  to  the  over-all  energy  budget  is  a  net  loss  of  energy 
by  the  atmosphere . 


V  THE  SHORT-WAVE  SOLAR  RADIATION  REGIME 


A .  GENERAL 

The  main  processes  effecting  the  distribution  of  solar  radiant 
energy  in  the  atmosphere  are  (l)  scattering  by  the  molecules  of  the  air, 
(2)  reflection  from  ground  or  clouds,  (s)  absorption  by  atmospheric 
constituents  such  as  ozone  and  water  vapor,  and  (4)  scattering  and 
absorption  by  particles  such  as  haze  and  dust.  For  the  present  case 
study,  the  first  two  of  these  were  studied  by  the  use  of  the  exact 
theory  of  radiative  transfer,  as  originally  derived  by  Chandrasekhar® 
and  extensively  applied  by  Sekera  and  collaborators.*®  The  latest 
information  on  the  distribution  of  ozone  was  used  to  determine  the  rate 
of  heating  by  ultraviolet  absorption  by  ozone,  and  an  empirical  method 
derived  by  Korb,  et  al.,**  was  applied  for  determining  water-vapor 
absorption.  The  effects  of  haze  and  dust  were  evaluated  by  means  of  a 
three-layered  atmospheric  model  adapted  from  Waldram's*®  measurements. 

Quantitative  values  of  the  total  contribution  and  the  diurnal 
variation  of  each  process  were  computed  for  the  selected  synoptic  situ¬ 
ation  at  Omaha.  The  methods  and  results  are  discussed  in  the  following 
sections . 

B.  MOLECULAR  SCATTERING  AND  SURFACE  REFLECTIONS 

1 .  Introduction 

Solar  radiation  is  returned  to  space  by  the  earth-atmosphere  system 
principally  through  four  interdependent  processes,  namely  (l)  scattering 
by  gaseous  molecules  of  the  air,  (2)  reflection  from  the  underlying  land 
or  water  surface,  (3)  scattering  by  liquid  water  droplets  composing  clouds, 
and  (4)  scattering  by  aerosol  particles  such  as  haze  and  dust.  The 
contributions  of  the  first  three  of  these  processes  are  considered  quanti¬ 
tatively  in  this  section.  First  the  combined  effects  of  molecular 
scattering  and  surface  reflection  for  a  cloudless  and  dust-free  atmosphere 
are  evaluated  for  an  idealized  atmospheric  model.  Then,  the  reflectance 
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of  a  cloud-plus-Rayleigh-atmosphere  system,  In  which  the  top  of  a 
uniform  cloud  deck  of  reflectance  A  =  0.50  is  located  at  various 
altitudes,  is  similarly  examined. 

The  model  atmosphere  is  one  of  finite  optical  thickness  in  the 
vertical,  homogeneous  and  of  infinite  extent  in  the  horizontal,  and 
non-absorbing.  The  atmosphere  is  bounded  on  top  by  a  perfect  vacuum  and 
at  the  bottom  by  a  diffusely  reflecting  surface  of  arbitrarily  set 
reflectivity  A  such  that  0  <  A  <  1.  For  the  case  of  clouds,  the  lower 
boundary  is  the  top  of  the  cloud  deck;  otherwise  it  is  assumed  to  be  a 
diffusely  reflecting  surface  at  sea  level. 

2 .  Method  of  Calculation 

The  exact  solution  of  the  general  equation  of  radiative  transfer, 
as  obtained  by  Chandrasekhar,®  was  used  for  the  calculations.  The 
atmospheric  model  is  characterized  by  the  normal  optical  thickness 
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where  v  is  the  frequency  of  the  radiation,  is  the  frequency-dependent 
volume-scattering  coefficient,  and  z  is  the  vertical  coordinate.  We 
assume  ttFq  units  of  energy  per  unit  frequency  interval  and  unit  time  are 
incident  on  a  unit  area  normal  to  the  direction  of  propagation  of  the 
parallel  radiation  at  the  top  of  the  atmosphere,  the  direction  being 
specified  by  the  solar  zenith  angle  6^  and  azimuth  (|)q.  The  intensity, 

I,  of  the  scattered  radiation  emerging  from  the  top  of  the  atmosphere  is 
a  function  of  the  nadir  angle  0  and  azimuth  (|).  Then  the  total  upward 
flux  Dy  through  a  unit  horizontal  area  per  unit  frequency  interval  and 
unit  time,  is  given  by  an  integration  of  intensity  over  the  downward 
hemisphere: 
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If  we  assume  that  the  underlying  surface,  be  it  soil,  vegetation, 
or  other  material,  reflects  isotropically  a  fixed  fraction.  A,  of  the 
radiation  incident  on  it,  then  the  upward  stream  is  composed  of  both 
scattered  and  surface-reflected  radiation.  When  the  integration  of 
Eq.  (3)  is  carried  out  for  this  case  (see  Coulson'*),  we  obtain 
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Here  =  cos  ©q,  and  i  is  a  known  function  of  only  the  optical  thickness. 
The  functions  derived  originally  by  Chandrasekhar,®  are  func¬ 

tions  of  both  Hq  and  t.  Numerical  values  of  them  are  available.  For  this 
work  the  following  tabulations  "were  used: 

Sekera  and  Blanch:^*  t  =  1.0,  0.25,  0.15 

Sekera  emd  Ashburn:*®  t  =  0.10,  0.05,  0.02,  0.01,  0.003,  0.001 

Chandrasekhar  and  Elbert:*®  t  =  0.50 

Sekera  and  Ashburn  reported  some  computational  difficulties  for  t  =  0.003 
and  T  =  0.001,  and  there  is  some  question  about  the  accuracy  of  their 
tabulated  values.  The  comparisons  given  by  those  authors  indicate,  how¬ 
ever,  that  the  final  values  are  reliable  to  about  three  places,  which  is 
sufficient  for  the  present  computations,  and  the  results  obtained  here 
appear  to  be  internally  consistent  throughout.  Their  tabulations  were 
used,  therefore,  without  alteration. 

3 .  Effect  of  Surface  Reflectivity 

It  was  found  in  a  previous  investigation  (Coulson*®),  that  the 
reflectance  of  the  atmospheric  model  is  critically  dependent  on  the 
reflectivity  of  the  underlying  surface.  In  order  to  define  the  dependence 
more  explicitly,  the  monochromatic  reflectance,  given  by 

"v  =  -P- 


was  computed  for  eight  different  values  of  surface  reflectivity  for  a 
number  of  optical  thicknesses  and  solar  zenith  angles.  Curves  of  Ry  as 
a  function  of  surface  reflectivity  A  for  large  optical  thickness  are 


27 


given  in  Fig.  7,  the  parameter  being  solar  zenith  angle  0^.  The  dashed 
curve  represents  the  reflectance  the  system  would  have  if  the  atmosphere 
were  not  present.  It  is  seen  that  the  atmosphere  itself  scatters  back  a 

O 

large  portion  of  this  short-wavelength  radiation  (X  =  3120  A),*  particu¬ 
larly  at  the  larger  solar  zenith  angles.  A  value  of  Ry  =  0.75  is  found 
for  A  =  0  for  the  largest  ©qi  and  the  reflectance  of  the  system  is  every¬ 
where  greater  than  0.34.  Even  here,  however,  the  surface  reflected 
radiation  is  predominant  for  large  values  of  A.  For  A  =  1.0  all  of  the 
curves  must  converge  at  the  point  Ry  =  1.0,  since  we  assume  a  non- 
absoi’bing  atmosphere. 

The  equivalent  curves  for  moderate  optical  thickness  (t  =  O.IO)  and 
for  very  small  optical  thickness  (t  =  O.OOl)  are  shown  by  Figs.  8  and  9, 
corresponding  to  the  yellow  (546oX)  and  infrared  (17,000  A),  respectively. 

The  very  minor  role  played  by  atmospheric  scattering  compared  to  surface 
reflection  at  the  longer  wavelengths  is  evident .  This  wavelength  dependence 
is  a  feature  of  considerable  practical  importance  in  the  design  and  oper¬ 
ation  of  radiation-sensing  satellites. 

A  complete  tabulation  of  the  monochromatic  reflectance  Ry  of  the 
atmosphere-surface  system  for  ten  different  optical  thicknesses,  from 
which  those  of  Figs.  7,  8,  and  9  were  selected,  is  Included  as  an 
Appendix  of  this  report. 

4 .  Effect  of  Solar  Angle 

The  reflectance  of  the  surface-atmosphere  system  is  strongly  affected 
by  the  angle  of  incidence  Gq  of  the  incoming  radiation.  Curves  of  system 
reflectance  as  a  function  of  iJq  =  cos  6q  for  four  different  values  of 
surface  reflectivity  are  reproduced  in  Fig.  10.  The  data  were  obtained 
by  an  integration  of  Eq .  (5)  with  wavelength,  the  incident  energy  being 
assumed  to  have  the  same  spectral  distribution  as  that  of  extra-terrestrial  - 
sunlight.  The  spectral  range  covered  is  3100  1  <  X  <  25,000  A,  representing 


*  The  conversions  of  optical  thickness  to  wavelength  are  for  z  =  0  (sea 
level)  unless  otherwise  specified. 


28 


MONOCHROMATIC  ALBEDO  OF  PLANE  PARALLEL  ATMOSPHERE 


FIG.  7 

MONOCHROMATIC  REFLECTANCE  OF  MODEL  ATMOSPHERE-PLUS-SURFACE  SYSTEM 
AS  A  FUNCTION  OF  SURFACE  REFLECTIVITY  FOR  OPTICAL  THICKNESS  r=  1.00 
FOR  FIVE  DIFFERENT  SOLAR  ZENITH  ANGLES 
T  =  1.00,  corresponding  to  wavelength  X=  3120A  at  sea  level. 
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MONOCHROMA 


ALBEDO  FOR  NO  ATMOSPHERE 


FIG.  8 

MONOCHROMATIC  REFLECTANCE  OF  MODEL  ATMOSPHERE-PLUS-SURFACE  SYSTEM 
AS  A  FUNCTION  OF  SURFACE  REFLECTIVITY  FOR  OPTICAL  THICKNESS  t=  0.10 
FOR  FIVE  DIFFERENT  SOLAR  ZENITH  ANGLES  „ 

T'0.10,  corresponding  to  wavelength  X=  5460A. 


FIG.  10 

REFLECTANCE  OF  MODEL  ATMOSPHERE-PLUS-SURFACE  SYSTtM  AS  A  FUNCpON  OF  COSINE 
OF  SOLAR  ZENITH  ANGLE  FOR  SPECTRAL  INTERVAL  3100A<  X  <  25,OOOA  FOR  DIFFERENT 

VALUES  OF  SURFACE  REFLECTIVITY 


about  95^  of  the  incident  solar  energy,  and  integration  was  by  the 
trapezoidal  rule.  The  strong  dependence  of  system  reflectance  on  solar 
position  at  small  values  of  poi  particularly  pronounced  for  small  values 
of  A,  results  from  backscatter,  both  primary  and  multiple,  in  the  long 
atmospheric  pathlengths  existing  for  large  solar  zenith  distances.  The 
amount  of  backscatter  decreases  rapidly  with  increasing  Hq,  however, 
resulting  in  small  changes  of  system  reflectance  for  greater  than  about 
0.40.  The  dominating  influence  of  surface  reflectivity  is  evident  in  the 
curves . 
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To  show  the  effects  of  the  different  spectral  regions,  the  solar 
spectral  range  was  broken  into  the  ultraviolet  (3100  ^  <  X  <  4000  A), 
visible  (4000  X  <  X  <  7000  X),  and  infrared  (7000  X  <  X  <  25,000  X),  and 
reflectance  curves  were  determined  for  each  region  as  a  function  of  solar 
angle.  The  results  are  shown  in  Fig.  11  for  four  values  of  A.  It  is 
seen  that  the  reflectance  characteristics  change  significantly  from  the 
ultraviolet  to  the  infrared.  A  relatively  strong  dependence  on  solar 


FIG.  11 

REFlfCTANCE  OF  MODEL  ATMOSPHERE-PLUS-SURFACE  SYSTEM  AS  A  FUNCTION  OF  COSINE 
OF  SOUyi  ZENITH  ANgl£  WITH  SPECTRUM  DIVIDED  INTO  ULTRAVIOLET 
{3100A  <  X  <  4000A),  VISIBlf  (40(X)X  <  X  <  7000X),  AND  INFRARED 
(7000X  <  X<  25,OOOX) 
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angle  continues  throughout  for  the  ultraviolet,  the  one  exception  being 
for  A  =  0.80,  but  nowhere  do  the  slopes  show  extreme  values.  In  the 
infrared,  however,  the  rapid  decrease  of  reflectance  with  decreasing 
solar  zenith  angle  shows  that  Rayleigh  scattering  of  the  long  waves  is 
significant  only  for  very  long  atmospheric  pathlengths.  Fc  .•  A  =  0.80 
the  system  Is  a  good  reflector  for  all  spectral  regions. 

5.  Variation  with  Altitude  of  Reflecting  Surface 

Because  of  the  dependence  of  optical  thickness  on  altitude  as  well 
as  on  wavelength,  the  system  reflectivity  is  a  function  of  the  altitude 
Zg  of  the  reflecting  surface.  The  discussion  so  far  has  been  for  Zg  =0. 
Similar  computations  have  been  made  for  a  surface  of  reflectivity  A 
located  at  altitudes  of  Zs  =  2  km,  5  km,  10  km,  20  km,  and  oo .  The  last 
case  corresponds  to  the  surface  being  above  the  atmosphere,  so  that  the 
system  reflectance  is  identical  to  A.  Again,  the  incident  flux  and 
spectral  distribution  of  extra-terrestrial  sunlight  is  assumed,  and  the 

O  0 

range  covered  is  3100  A  <  X  <  25,000  A.  The  altitude-wavelength  dependence 
for  a  given  optical  thickness  was  taken  from  Deirmendjlan.  ^ Computations 
were  made  for  a  single  reflectivity,  A  =  0.50  having  been  chosen  as  a 
reasonable  value  for  the  reflectivity  of  an  actual  cloud. 

The  upward  flux  as  a  function  of  wavelength  for  Gq  =  66.4®  is  shown  in 
Fig.  12,  for  the  three  altitudes  Zg  =  0,  5  km,  oo  .  Ordinate  values  for 
the  curve  Zg  =  oo  are  simply  one-half  those  of  the  incident  radiation, 
since  the  atmospheric  effect  is  nil  for  this  case.  Comparison  of  the 
curves  shows  that  the  effect  of  atmospheric  scattering  by  the  air  above 
a  highly  reflecting  cloud  top  is  minor,  or  completely  negligible,  at  this 
sun  elevation  for  the  entire  region  X  >  8000  In  the  visible  and  ultra¬ 
violet,  however,  the  dependence  of  upward  flux  on  Zg  is  important,  the 
value  being  41^  higher  for  Zg  =  0  than  for  Zg  =  oo  at  X  =  3100  X.  The 
effect  is  less  at  smaller  solar  zenith  angles,  but  the  corresponding 
increment  is  still  13^  for  0©  =  0. 

The  case  for  0q  =  88.9®  is  given  by  Fig.  13,  in  which  curves  for 
five  different  values  of  Zg  are  included.  A  major  change  with  Zg  is  seen 
here,  and  the  influence  extends  well  into  the  infrared.  Even  for 
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FIG.  12 

UPWARD  FLUX  FROM  THE  TOP 
OF  MODEL  ATMOSPHERE  AS  A 
FUNCTION  OF  WAVELENGTH  FOR 
SOLAR  ZENITH  ANGLE  66.4°  WITH 
A  SURFACE  OF  REFLECTIVITY  0.50 
LOCATED  AT  DIFFERENT  ALTITUDES 


Zs  =  20  km,  atmospheric  scattering  increases  the  upward  flux  by  46^  over 
its  value  for  Zg  =  oo  at  X  =  3100  A. 

If  the  spectral  range  is  again  divided  into  ultraviolet,  visible, 
and  infrared  regions  and  the  angular  integrations  described  above  are 
performed,  it  is  a  simple  matter  to  compute  the  reflectance  of  the  system 
as  a  function  of  n©*  'jFh®  result  of  this  procedure  fop  the  ultraviolet, 
for  A  =  0.50,  is  shown  in  Fig.  14.  For  the  sun  near  the  horizon 
(Gq  =  88.9°) the  reflectance  of  the  system  changes  slowly  with  Zg ,  de¬ 
creasing  only  to  0.70  for  Zg  =  20  km.  The  region  of  maximum  sensitivity. 
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FIG.  13 


UPWARD  FLUX  FROM  THE  TOP  OF  MODEL 
ATMOSPHERE  AS  A  FUNCTION  OF  WAVELENGTH 
FOR  SOLAR  ZENITH  ANGLE  88.9*  WITH  A  SURFACE 
OF  REFLECTIVITY  0.50  LOCATED  AT  DIFFERENT 
ALTITUDES 
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FIG.  14 

SYSTEM  REFLECTANCE  AS  A  FUNCTION  OF  ALTITUDE  OF  SURFACE 
OF  REFLECTIVITY  0.50,  FOR  THE  ULTRAVIOLET (3 lOoA  <  X<  4000A) 

AT  FOUR  DIFFERENT  SOLAR  ZENITH  ANGLES 

indicated  by  the  greatest  slope  of  the  curves,  decreases  in  altitude  with 
decreasing  solar  angle.  For  small  zenith  angles  the  sensitivity  is  small 
at  all  altitudes  because  of  the  dominance  of  surface-reflected  radiation. 

The  equivalent  curves  for  the  visible  and  infrared  are  reproduced  in 
Fig.  15.  The  dominant  role  in  surface-reflected  radiation  makes  the 
reflectance  insensitive  to  Zg  for  the  infrared,  but  for  the  visible  a 
large  dependence  is  indicated  at  low  sun  elevations. 
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FIG.  15 

SYSTEM  REFLECTANCE  AS  A  FUNCTION  OF  ALTITUDE  OF  SURFACE  OF 
REFLECTJVITY  0.50,  VISIBLE  (4000A  <  X  <  7000A)  AND  INFRARED 
t7000A  <  X<  25,000A)  AT  FOUR  DIFFERENT  SOLAR  ZENITH  ANGLES 


6.  Discussion 


From  a  practical  standpoint,  it  appears  that  the  wavelength  dependence 
of  the  light  emerging  from  the  top  of  the  atmosphere  could  be  used  as  a 
measure  of  the  amount  of  Rayleigh  scattering  occurring  above  a  cloud  top, 
or,  equivalently,  the  altitude  of  the  cloud  top.  In  order  to  assess  the 
possibilities,  we  should  take  finother  look  at  the  assumptions  underlying 
these  computations. 
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First,  the  reflectance  of  the  cloud  top  is  assumed  to  be  0.50, 
independent  of  wavelength.  This  value  of  A  was  arbitrarily  chosen  as 
being  reasonable,  and  not  as  realistic  for  any  given  cloud  type,  but  the 
character  of  the  curves  would  persist  through  a  reasonable  change  of  A. 
Furthermore,  wavelength  independence  of  A  is  not  required,  but  only  that 
the  functional  dependence  be  known.  From  theory,  one  would  expect  the 
wavelength  dependence  of  cloud-droplet  scattering  to  be  minimal,  since 
the  cloud  droplets  are  well  into  the  Mie  scattering  domain  and  consist 
of  a  distribution  of  sizes.  This  is  qualitatively  borne  out  for  the 
visible  spectrum  by  the  lack  of  color  in  cloud-reflected  or  transmitted 
sunlight . 

The  second  assumption,  that  of  a  horizontally  stratified  atmosphere 
resting  on  a  flat  earth.  Introduces  an  artificiality  that  must  be  removed 
for  hemispheric  flux  measurements  from  a  satellite,  to  be  useful  for 
cloud  height  determination.  Approximation  of  atmospheric  sphericity  by 
many  plane  areas  of  different  orientations  will  probably  be  satisfactory, 
however,  for  present  purposes. 

Third,  radiative  scattering  by  dust  and  haze  has  not  been  considered 
here.  A  first  approximation  to  the  effect  could  be  obtained  by  the  method 
suggested  by  Sekera^®  if  the  physical  properties  of  the  particles  and 
their  spatial  distribution  were  known.  These  quantities,  however,  are 
extremely  variable  in  the  atmosphere,  and  data  are  sparse.  Pending 
further  work  in  this  area,  the  effects  of  such  atmospheric  turbidity  must 
serve  only  to  Increase  the  uncertainties  Inherent  in  the  method. 

It  would  appear  from  such  preliminary  considerations  that  this  method 
of  cloud-height  determination  is  sufficiently  attractive  to  merit  more 
detailed  developmental  work. 

Two  other  possibilities  for  cloud-height  determination  by  light- 
scattering  measurements  should,  perhaps,  be  mentioned;  (l)  Instead  of 
integrating  the  intensity  over  the  hemispheric  solid  angle,  as  suggested 
above,  there  are  some  important  advantages  in  measuring  the  intensity  of 
the  light  scattered  into  a  small  solid  angle  in  a  given  direction. 

(2)  The  extremely  different  polarizing  properties  of  air  molecules  and 
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cloud  particles  might  be  used  to  determine  the  optical  thickness  of  the 
air  overlying  a  cloud  deck,  and  hence  the  height  of  the  cloud  top. 

7 .  Application  to  Synoptic  Situation 

The  above  results  were  applied  to  the  particular  synoptic  situation 
under  investigation,  yielding  values  of  the  upward  flux  of  solar  energy 
from  the  top  of  the  atmosphere  and  the  energy  absorbed  by  the  atmosphere- 
ground  system  at  Omaha,  as  a  function  of  time.  Clouds  were  assumed  to 
reflect  0.50  of  the  radiation  incident  on  them,  and  for  partly  cloudy  sky 
the  data  were  weighted  by  the  fractional  sky  cover.  The  reflectivity  of 
the  snow-covered  ground  was  taken  as  A  =  0.80,  and  for  bare  ground 

O 

A  =  0.10.  No  correction  was  made  for  the  wavelength  regions  X  >  25,000  A 
and  X  <  3100  A.  Since  those  combined  regions  contain  only  about  5^  of  the 
solar  energy  and  constitute  both  strongly  scattered  and  weakly  scattered 
portions,  their  effects  on  the  curves  would  be  small. 

The  solid  curves  of  Fig.  16  show  the  flux  of  solar  energy  returned 
to  space  by  the  atmosphere-ground  system  at  Omaha  for  the  four  days 
March  12-15,  1957,  while  the  dashed  curves  give  the  energy  absorbed  by 
the  ground  during  the  same  period.  The  great  effedt  of  surface  re-‘ 
fleet! vlty  is  evident,  being  particularly  pronounced  for  the  case  of  snow 
cover  during  the  forenoon  of  March  15.  Cloud  height  has  a  minor  effect 
as  long  as  the  cloud  reflectivity  is  independent  of  height.  This  can  be 
seen  by  comparing  the  curves  of  March  13  (cloud  top,  8300  m)  with  those 
of  March  14  (cloud  top,  3000  m  and  5700  n) . 

C.  ABSORPTION  BY  WATER  VAPOR 

Hie  absorption  of  solar  energy  by  the  usual  amounts  of  atmospheric 
water  vapor  are  not  unimportant  in  the  atmospheric  heat  balance.  Korb, 
et  £l.,^^  found  the  amount  of  absorbed  energy  to  vary  from  178  ly/day  for 
high  water  vapor  content  and  high  incident  energy  to  zero  for  the  polar 
night.  Values  of  5  to  10^  of  the  incident  energy  were  common.  Kimball^® 
gave  values  of  6  to  125^  for  moderate  water-vapor  pathlengths . 

By  processing  the  absorption  data  of  Howard,  Burch,  and  Williams,®® 
Korb"  found  that  the  amount  of  solar  energy  Ey  absorbed  by  water  vapor  in 
the  entire  thickness  of  the  atmosphere  can  be  approximated  by 
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UPWARD  FLUX  OF  SOLAR  RADIATION  FROM  THE  TOP  OF  THE 
AS  A  FUNCTION  OF  LOCAL  TIME  FOR  A 


ENERGY  ABSORBED  BY  THE  GROUND 


OMAHA  LOCAL  TIME  OMAHA  LOCAL  TIME 

MARCH  13  MARCH  14 


FIG.  16 

FLUX  OF  SOLAR  RADIATION  FROM  THE  TOP  OF  THE  ATMOSPHERE  AND  SOLAR  ENERGY  ABSORBED  BY  GROUND  AT  OMAHA 
AS  A  FUNCTION  OF  LOCAL  TIME  FOR  A  SELECTED  SYNOPTIC  SITUATION  IN  MARCH  1957 
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FIG.  16 

OF  THE  ATMOSPHERE  AND  SOLAR  ENERGY  ABSORBED  BY  GROUND  AT  OMAHA 
FOR  A  SELECTED  SYNOPTIC  SITUATION  IN  MARCH  1957 
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log  Ey  =  -0.740  +  0.347  log  V  -0.056  (log  v)**  -0.006  (log  v)®  (6) 

where  V  Is  the  preclpltable  water  (cnj)  in  the  path  of  the  radiation.  If 
V  is  the  precipitable  water  in  the  vertical  column,  then  for  the  incident 
parallel  solar  radiation  Vj  =  v  sec  Qq,  and  for  diffuse  radiation, 

Vj)  1 .66  V . 

1.  Absorption  for  Clear-Sky  Conditions 

If  for  clear  skies,  is  the  energy  absorbed  from  a  vertically 
incident  beam  per  unit  time  by  the  water  vapor  in  a  vertical  air  column 
of  a  unit  cross  section,  then  we  can  write 

®vC  “  ^IC  ®SC  ^DC 


where  Ejq  is  the  energy  absorbed  from  the  incident  parallel  sunlight,^ 

EgQ  is  that  absorbed  from  the  scattered  sunlight,  and  Ej^  is  the  energy 
absorbed  from  the  light  diffusely  reflected  from  the  underlying  surface. 

The  term  Ej^  generally  represents  the  major  component  of  E^q  and  is 
easily  computed  by  Eq .  (6)  from  the  data  of  v  and  0q.  For  practical  pur¬ 
poses,  EgQ  is  negligible  in  comparison  to  the  other  terms  since  water 
vapor  absorbs  .mainly  in  the  Infrared  region  for  which  molecular  scattering 
is  very  small.  The  term  Ej)(j  is  more  important,  especially  for  highly 
reflecting  surfaces.  The  total  water  vapor  pathlength  encountered  by  the 
reflected  radiation  is  Vr  =  Vi  +  Vp,  the  two  terms  corresponding  to  the 
downward  and  upward  traverses  of  the  atmosphere,  respectively.  But  the 
absorption  in  the  downward  traverse  was  already  accounted  for  in  Ejq. 

Then,  for  a  surface  of  reflectivity  A,  we  have  for  the  diffusely  reflected 
radiation  in  a  clear  atmosphere 

Edc  =  A[E(V„)  -  Ejc]  (8) 


Where  e(Vr)  is  evaluated  by  Eq.  (6)  for  the  pathlength  Vr.  If  the  change 
of  incident  energy  with  solar  angle  ©q  is  taken  into  account,  the  total 
energy  flux  absorbed  by  water  vapor  in  the  clear  atmosphere  is 


=  E 


vC 


cos  ©o  =  [Ejc  +  A[e(Vr)  -  Ejcj}  cos  ©^  . 


(9) 
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2.  Absorption  for  Overcast  Conditions 


The  analysis  for  overcast  conditions  is  complicated  by  the  presence 
of  clouds.  A  fraction  A^i  of  the  Incident  radiation  is  reflected  (actually 
back-scattered)  from  the  cloud  top,  and  therefore  undergoes  two  traverses 
of  the  water  vapor  above  the  cloud.  By  reasoning  similar  to  that  leading 
to  Eq.  (9),  it  is  seen  that  the  amount  of  energy  absorbed  above  the  cloud 
top  is 

Ea  =  §Ia  +  At  [E(VRa)  -  Eiaj}  .  (lO) 

Here  VRa  is  the  total  pathlength  for  the  reflected  radiation  above  the 
cloud  top,  the  downward  traverse  being  for  parallel  and  the  upward 
traverse  for  diffuse  radiation.  A  small  error  would  be  introduced  by 
the  actual  pathlength  being  different  from  since  some  of  the  radiation 
penetrates  a  distance  into  the  cloud  before  being  scattered  back,  and  some 
multiple  scattering  takes  place  Inside  the  cloud.  Neiburger's  measurements 
indicate  that  the  relatively  thin  clouds  are  efficient  backscatterers,  a 
cloud  of  reflectance  Aj  =  0.50  being  generally  less  than  500  feet  thick. 
Since  much  of  the  back-scattered  radiation  would  never  reach  500  feet  into 
the  cloud,  300  feet  is  perhaps  a  reasonable  average  value.  The  cloud  top 
height  would  thus  be  effectively  decreased  by  this  amount. 

A  small  part  of  the  fraction  1  -  A.ji  which  is  transmitted  through  the 
cloud  top  i'  ■'bsorbed  by  the  liquid  water  of  the  cloud  droplets.  The 
effect  of  this  absorption  on  the  water  vapor  problem  is  second-order, 
however,  anu  will  be  neglected  here. 

The  radiation  entering  the  cloud  quickly  becomes  diffuse,  and  the 
effective  pathlength  for  the  region  below  the  cloud  top  is  changed  thereby 
from  V]3  sec  0^  to  1.66  Vj^,  where  Vj^  is  the  precipitable  water  in  vapor 
form  below  the  cloud  top.  Multiple  scattering  by  the  droplets  inside  the 
cloud  would  increase  the  mean  pathlength  somewhat  further.  Korb,  et  al . , ^ ^ 
derived  a  "detour  factor"  to  account  for  this  path  extension,  the  value  of 
which  indicates  that  there  may  be  up  to  a  50^  increase  for  certain  con¬ 
ditions.  Smallei  increases  are  more  common,  however.  Their  result  of  the 
detour  factor  being  less  than  unity  seems  of  questionable  significance. 
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and  they  point  out  that  their  results  represent  only  order-of-'inagnitude 
corrections.  Until  more  advanced  methods  are  available  and  more  infor¬ 
mation  on  the  cloud  characteristics  themselves  is  obtained,  it  seems 
premature  to  Introduce  a  refinement  for  this  pathlength  Increase  by 
multiple  scattering  inside  the  cloud. 

For  the  present  computations,  it  is  assumed  that  the  cloud  effect 
for  water-vapor  absorption,  in  the  region  below  thev  cloud  top,  is  a 
decrease  of  the  downward  flux  due  to  the  reflectivity  A.j.  of  the  cloud 
top,  and  a  change  of  the  effective  pathlength  by  a  factor  of  1.66  by  the 
radiation  becoming  diffuse.  The  total  pathlength  of  the  radiation  trans¬ 
mitted  to  the  surface  is  Vj  =  v^^  sec  +  1.66  v^^,  where  v^  is  the  pre- 
clpl table  water  vapor  above  the  cloud  top  and  v^^  is  that  below.  The 
fraction  of  the  incident  energy  traversing  this  distance  is  approximately 
1  -  Ax.  It  is  easily  seen  that  the  energy  absorbed  below  the  cloud  top 
from  the  downward  stream  is 

=  (1  -  At)[e(Vj)  -  E(Va)]  .  (11) 

There  will  also  be  some  absorption  of  the  surface-reflected  radiation. 

The  additional  pathlength  for  this  diffuse  component  is  V  -ss!  1.66  v  and 
the  absorbed  energy  is 


Ej.  =  At[e(Vj.)  -  E(Vi)]  .  (12) 

The  total  energy  absorbed  by  water  vapor  for  the  overcast  case  for 
solar  zenith  angle  0^  is  then 

^0  =  So  .  (13) 

The  rate  of  energy  absorption  by  the  atmospheric  water  vapor  at 
Omaha  was  computed  by  Eq.  (6)  for  the  synoptic  case  of  March  12-15,  1957. 
Typical  results  are  shown  as  a  function  of  local  time  in  the  diagrams  of 
Fig.  17.  The  actual  radiosonde  humidity  measurements  were  used  to  compute 
the  various  pathlengths  when  such  data  were  available,  and  at  other  times 
humidity  data  from  the  vertical  time  section  were  used.  For  periods  of 
partial  cloudiness  the  computations  were  carried  out  for  both  clear  and 
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FIG.  17 

RATI  OF  SOLAR  ENERGY  ABSORPTION  BY  ATMOSPHERIC  WATER  VAPOR 
AT  OMAHA  AS  A  FUNCTION  OF  LOCAL  TIME  FOR  SELECTED 
SYNOPTIC  SITUATION 
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overcast  conditions  and  the  final  results  were  weighted  according  to  the 
cloud  amount.  The  discontinuities  of  the  curves  are  the  result  of  assumed 
discontlnuously  changed  conditions  of  water-vapor  content,  cloudiness,  or 
surface  reflectivity.  The  cloud  reflectivity  was  taken  as  =  0.50  for 
overcast.  Surface  reflectivity  was  again  assumed  to  be  A  =  0.10  for  the 
snow- free  surface  and  A  =  0.80  for  new  snow. 


For  clear  conditions  with  snow-free  ground  (March  15,  P.M.)  practi¬ 
cally  all  of  the  absorbed  energy  is  from  the  downward  traverse  of  the 
atmosphere,  the  maximum  rate  of  absorption  being  of  the  order  0.13  -  0.15 
ly/min.  Absorption  of  surface-reflected  radiation  contributes  only  about 
1^  or  2^  for  the  snow- free  case,  but  strongly  reflecting  new  snow 
(March  15,  A.M.)  increases  the  surface-reflected  contribution  to  a  maximum 
of  about  20^  of  the  total.  The  somewhat  greater  total  absorption  which 
occurred  on  March  14  is  the  result  of  a  higher  total  water-vapor  content. 

A  change  of  cloud-top  height  from  2.7  km  to  5.5  km  explains  the  shift  of 
relative  magnitudes  of  the  dashed  curves  at  about  0900  on  March  14.  The 
cloud-reflected  radiation  which  is  absorbed  above  the  cloud  top,  repre¬ 
sented  by  the  bottom  curve  of  the  March  14  diagram,  reaches  a  maximum  of 
about  6^  of  the  total. 


Integration  of  the  time  rate  of  energy  absorption  from  sunrise  to 
sunset  yields  the  total  energy  absorption  and  percent  of  extraterrestrial 
incident  flux  for  each  day  of  the  period  March  12-15  given  in  Table  I. 

It  is  seen  that  the  values  agree  well  with  those  obtained  by  other 
investigators . 


TABLE  I 


TOTAL  SOLAR  ENERGY  AND  PERCENT  OF  EXTRATERRESTRIAL 
INCIDENT  FLUX  ABSORBED  BY  ATMOSPHERIC  WATER  VAPOR 
AT  OMAHA  DURING  EACH  DAY  OF  THE  PERIOD  MARCH  12-15,  1957 


Date 

Total  Absorbed 
Solar  Energy 
(langleys) 

Percent  of  Extraterrestrial 
Incident  Flux 

March  12 

66.3 

10.8 

March  13 

65.2 

10.6 

March  14 

77.5 

12.6 

March  15 

56.9 

9.3 
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McDonald^*  used  the  well-known  absorption  data  of  Fowle  to  obtain 
an  empirical  absorptivity  a  of  the  form  a(u)  =  0.077(u)®*®°  where  u(cm) 
is  the  total  preclpitable  water  vapor  in  the  optical  path.  Application 
of  this  method  to  the  data  of  March  12  for  Omaha  yields  an  absorption 
value  of  roughly  20^^  less  than  that  of  Table  I.  McDonald  points  out  that 
a  question  has  been  raised  as  to  the  validity  of  Fowle 's  data,  and  con¬ 
cludes  that  absolute  accuracies  of  atmospheric  water-vapor  absorption 
calculations  of  not  better  than  30^  are  possible  at  present.  In  view, 
however,  of  the  independent  absorption  measurements  of  Howard,  et  al.,^® 
which  were  used  byKorb,^^  the  30^  maximum  accuracy  appears  pessimistic, 
but  the  errors  have  not  been  well  evaluated  as  yet  in  the  present  study. 

D.  EFFECTS  OF  HAZE  AND  DUST 

The  radiative  effects  of  haze,  dust,  and  other  such  particulates  in 
the  atmosphere  are  not  well  known,  partly  because  of  lack  of  information 
on  the  particles  themselves  and  partly  because  of  the  difficulties  in¬ 
herent  In  solving  the  radiative  transfer  problem  in  a  turbid  atmosphere. 
Measurements  of  the  physical  characteristics  of  atmospheric  suspensolds 
emphasize  their  extreme  variability  of  size-frequency  distributions  and 
concentrations,  both  of  which  are  time-  and  space-dependent.  A  simple 
power  law  is  commonly  used  as  the  best  description  of  the  average  size- 
frequency  distribution,  while  Deirmendjian^*  combined  similar  laws  with 
an  exponential  decrease  of  the  concentration  with  height  for  computing 
the  intensity  distribution  in  the  solar  aureole.  The  measurements  further 
indicate  characteristic  differences  between  the  aerosols  found  in  conti¬ 
nental  as  opposed  to  maritime  air  masses,  and  urban  pollutants  are  different 
from  those  of  predominantly  rural  areas. 

The  problem  of  radiative  transfer  in  a  turbid  medium  is  receiving 
considerable  emphasis  by  various  authors,  but  no  completely  satisfactory 
solution  for  the  general  problem  is  available  as  yet.  The  exact  solution 
of  the  transfer  problem  for  molecular  atmospheres,  as  obtained  originally 
by  Chandrasekhar®  and  further  extended  by  Sekera  and  collaborators,'®*'® 
is  a  major  recent  development  in  the  field,  and  it  has  yielded  quantitative 
data  describing  completely  the  radiation  field  at  the  upper  and  lower 
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boundaries  of  realistic  atmospheric  models  (Coulson,  et  al.^®).  Sekera*® 
has  shown  a  way  of  introducing  large-particle  scattering  into  the 
appropriate  equations  to  give  a  first  approximation  for  the  effects  of 
turbidity,  and  Deirmendjian**  has  used  the  method  for  computing  the 
radiative  distribution  near  the  forward  scattering  angles.  Fraser®^  has 
used  the  Mie-coef ficlents  to  compute  the  normalized  scattering  matrix  for 
three  different  models  of  atmospheric  aerosols,  thereby  providing  the 
possibility  of  solving  the  radiative  transfer  equation  for  a  turbid 
atmosphere.  The  numerical  solution,  however,  is  still  Intractable,  even 
for  modern  computing  equipment,  because  of  the  complicated  dependence  of 
aerosol  scattering  on  scattering  angle.  The  basic  Mle  theory  provides 
for  absorption  also,  but  further  computational  difficulties  arise  from 
its  Introduction  into  the  system, 

For  the  present  problem  of  determining  the  radiant  energy  budget  of 
the  selected  synoptic  situation  the  measurements  by  various  authors  may 
be  used.  Most  of  the  work  has  boon  for  the  visible  region  of  the  spectrum, 
but  the  approximate  independence  of  aerosol  scattering  with  wavelength 
makes  it  feasible  to  use  the  measurements  for  the  whole  solar  energy 
spectrum  for  a  rough  evaluation  of  the  energy  budget. 

Waldram^®  determined  the  scattering  coefficient  Og  and  absorption 
coefficient  for  the  visible  in  both  clean  and  Industrially  polluted 
atmospheres.  The  scattering  coefficient,  given  by 


a 


s 


i  ^ 

F  dx 


(14) 


where  F  is  flux  and  x  is  distance,  was  obtained  by  integrating  the 
intensity  Ig  of  the  scattered  light  over  the  entire  solid  angle  Q.  Thus 
the  energy  attenuated  from  the  original  incident  flux  was 


-  dF  -  a,  r  dx  ■  f  Ig  du  .  (15) 

Values  of  Og  varied  from  8.0  x  for  very  clean  air  to  2.5  x  10“®/ft 

for  foggy  conditions  of  1000  yards  visual  range.  In  moderate  visibilities 
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(about  6  miles)  Og  si  1.6  x  10“^/ft.  Considerable  variation  with  height 
was  found  to  exist  through  the  entire  layer  0  to  30,000  feet  altitude. 

In  clean  atmospheres  there  was  generally,  but  not  always,  a  decrease  of 
as  with  height  of  an  amount  roughly  proportional  to  the  decrease  of  pure 
air  scattering  with  height.  For  polluted  atmospheres  large  fluctuations 
occurred,  particularly  in  the  lower  atmospheric  layers. 

The  shape  of  the  polar  scatter  diagram  is  of  obvious  importance  in 
energy  budget  studies.  For  primary  molecular  scattering,  the  scattered 
energy  is  directed  equally  into  the  upward  and  downward  hemispheres, 
regardless  of  sun  elevation.  Multiple  molecular  scattering  changes  the 
partition  a  small  amount  by  increasing  the  backscatter  to  somewhat  more 
than  half  the  total  (Coulson^?).  The  predominantly  forward  scattering 
of  aerosol  particles,  however,  alters  this  simple  pattern,  making  the 
downward  component  larger  in  general,  and  introducing  a  solar  angle 
dependence.  Only  for  the  sun  at  the  horizon  would  the  upward-  and 
downward-scattered  components  be  equal.  Various  empirically  derived 
curves  have  been  fitted  to  the  measured  polar  intensity  distributions, 
but  they  are  of  dubious  value  because  of  the  great  variability  of  the 
natural  aerosol  properties  and  the  difficulty  of  measuring  accurately 
the  intensity  near  the  forward  direction.  The  more  fundamental  approach, 
employed  by  Fraser,*^  of  comparing  the  measured  polar  distribution  with 
that  computed  for  representative  aerosol  models  of  known  characteristics 
holds  great  promise,  but  further  work  is  required.  For  instance,  Fraser 
found  wide  discrepancies  between  his  computed  distributions  and  Waldram's 
measurements.  He  suggested  an  aerosol  model  with  a  larger  number  of 
particles  with  high  refractive  index,  but  apparently  such  a  model  has  not 
been  investigated. 

From  measurements  of  total  transmission,  T,  and  the  scatter  coef f iclent ,  as, 
it  is  possible  to  compute  the  absorption  coefficient,  a^,  from  the  relation 

agX  -  In  T  (16) 


where  x  is  distance  through  the  medium.  Waldram  concluded  from  his 
measurements  that  the  "phenomena  of  absorption  and  scatter  appear  to  be 
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independent.”  It  Is  known  from  theory  that,  for  a  single  particle  at 
least,  the  two  are  not  independent,  but  apparently  the  effects  of  inter¬ 
dependence  tended  to  be  small  for  the  conditions  prevailing  at  the  time 
of  Waldram's  measurements.  Generally  the  absorption  was  found  to  be 
appreciable  only  when  smoke  was  present,  at  which  time  was  usually  in 
the  range  1.0  x  10“*/ft  <  <  ^-0  x  10“*/ft  and  was  confined  to  more  or 

less  distinct  layers  of  a  few  hundred  feet  thickness.  In  smoke- free 
conditions  Oa  was  of  the  order  of  10”®/ft,  and  in  a  case  of  heavy  pollution 
it  reached  6  x  10”*/ft. 

For  a  rough  evaluation  of  the  effect  of  scattering  and  absorption 

in  the  March  situation  at  Omaha,  a  three-layer  model  of  the  atmosphere 

was  constructed  and  the  parameters  of  each  layer  were  estimated  from 
Waldram's  measurements.  The  model  is  shown  schematically  in  Fig.  18. 

The  bottom  layer,  2000  feet  in  thickness,  is  assumed  to  have  variable 
scattering  and  absorption  coefficients,  the  values  depending  on  the  amount 
of  industrial  pollution  present.  Low  visibilities  at  Omaha  on  March  13 
indicate  the  presence  of  industrial  pollution,  while  very  clear  air  on 
March  15,  following  the  frontal  passage,  shows  no  such  pollution.  The 
coefficients  were  assumed  to  change  accordingly,  the  assumed  values  being 
Cg  =  7.5  X  10“®/ft,  Og  =  1.0  X  10“*/ft  for  polluted  air  and 

Cg  =  1.0  X  10“®/ft,  cjg  =  0  for  Clear  air.  (The  effects  of  haze  and  dust 

are  not  computed  for  cloudy  periods,  as  their  effects  are  probably  small 
compared  to  those  of  clouds.)  The  assumption  of  negligible  absorption  in 
the  absence  of  industrial  pollution  is  substantiated  by  both  the  measure¬ 
ments  of  Waldram  and  an  equivalent  statement  by  Van  de  Hulst.^®  The  layer 
between  2000  and  30,000  feet  altitude  is  assumed  to  have  no  absorption 
and  an  average  scattering  coefficient,  derived  from  Waldram's  measurements 
for  clear  air,  of  Og  =  1  x  10“®/ft.  Above  30,000  feet  Rayleigh  scattering 
is  assumed,  an  optical  thickness  of  t  =  0.04  being  reasonable  for  the 
entire  mass  of  air  above  that  level.  Optical  thicknesses  due  to  scattering 
in  the  first  and  second  layers  were  computed  by  the  definition 

^3(2)  =  \  dz 

Jzi 
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FIG.  18 

SCHEMATIC  DIAGRAM  OF  THREE-LAYER  ATMOSPHERIC  MODEL 
Adapted  from  Waldram's  measurements 
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and  similarly  for  the  absorption  optical  thickness,  Finally,  the 

attenuation  of  the  incident  parallel  beam,  inclined  at  angle 
Gq  =  arc  cos  from  the  vertical,  by  scattering  in  the  entire 
atmosphere  is 

1  . 

For  convenience  of  computation  it  was  assumed  that  the  absorption  and 
scattering  are  independent  processes  and  that  the  absorption  occurs  from 
the  final  surface  radiation  field  existing  after  scattering  (i.e., 
scattering  occurs  first).  Thus  it  is  reasonable  to  consider  all  of  the 
radiation  as  being  diffuse  for  the  absorption  process,  and  that  the 
effective  absorption  optical  thickness  is  thereby  1.66  independent 
of  solar  angle.  It  is  readily  seen  that  the  following  relations  exist 
among  the  transmissivity  T,  absorptivity  a,  and  backscatter  b  for  the 
total  atmospheric  model: 

1  =  T  +  a  +  b 

T  =  e-i.66T&  [e’^^s/uo  +  f(i  -  e"'^s^o)] 

/  / 

a  =  (1  -  e"‘*ee'^a)[6"^s/|io  +  f(i  _  e"'^s/^Io)] 

b  i  (1  -  f)(l  -  . 

Here  f  is  the  fraction  of  the  scattered  energy  which  is  directed  into 
the  downward  hemisphere. 

To  estimate  the  portion  of  the  energy  which  is  backscattered  and 
lost  to  the  system  it  is  necessary  to  integrate  the  appropriate  intensity 
phase  function  over  the  upward  and  downward  hemispheres.  Plaskovskaia- 
Fesenkova*®  describes  the  scattering  diagram  for  clear  air  by  the  empiri¬ 
cally  derived  function 

P(\lr)  =  1  +  0.49  cos  \jr  +  1.093  cos^  ijr  +  0.704  cos®  \|a  (18) 

where  \|r  is  the  scattering  angle  measured  from  the  forward  direction.  The 
forward  component  of  the  scattered  energy  is,  except  for  a  constant, 
given  by 
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P(\lr)  cos  tIt  sin  ilr  dijf 


(19) 


t/s 


and  the  backward  component  by 


cos  ilr  sin  \lr  dilr 


(20) 


On  carrying  out  the  integrations  and  taking  the  ratio  of  the  components, 
we  find  the  forward  component  to  be  1.51  times  the  backward  component. 

Gibbons*'’’  describes  the  scattering  diagram  for  a  clear  atmosphere 
by  the  empirically  derived  function 

Gi(\lr)  =  1  +  coflis  .  (21) 


For  this  the  forward  component  is  1.36  times  the  backward  component. 

These  values  would  be  the  ratio  of  downward  to  upward  energy  in  the 
atmosphere  only  for  vertical  incidence.  Obviously  for  horizontal  inci¬ 
dence  the  ratio  would  be  1.0,  and  an  effective  value  applicable  for  the 
whole  day  would  be  somewhere  between  the  two  limiting  cases.  For  present 
purposes  a  value  of  1.2  is  taken  as  a  reasonable  ratio  for  the  clear 
atmosphere . 

For  a  noticeably  hazy  atmosphere  Gibbons*'^  fitted  the  observed 
intensity  distributions  adequately  by  the  relation 

GgCilf)  =  1+9  cos’e  .  (22) 

Substitutions  of  this  in  Eq.  (19)  and  Eq.  (20)  yields  a  ratio  of  forward- 
to-backward-scattered  energy  of  4.20,  applicable  for  vertical  incidence. 

On  considering  the  change  with  sun  elevation,  a  ratio  of  2.0  seems 
reasonable  for  the  effective  mean  value  for  the  polluted  atmosphere  at 
Omaha.  This  predominantly  forward  scattering,  however,  prevailed  in  only 
the  lower  2,000  feet,  which  contains  only  about  8^. of  the  total  atmosphere, 
so  the  mean  scattering  pattern  would  not  be  changed  greatly  from  the  clear 
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case.  A  ratio  of  forward-to-backward-scattered  energy  of  1.5  seems 
reasonable  for  the  total  atmosphere  with  the  presence  of  the  low-level 
pollution . 

The  fractional  distribution  of  the  extraterrestrial  solar  flux 
computed  for  the  above  model  is  shown  in  Fig.  19  as  a  function  of  local 
time,  for  the  clear  atmosphere  of  March  15  (solid  curves)  and  the  in¬ 
dustrially  polluted  atmosphere  of  March  13  (dashed  curves).  The  existence 
of  a  2,000-foot  pollution  layer  had  little  effect  on  the  total  back- 
scattering  characteristics.  The  effect  of  the  larger  scattering  co¬ 
efficient  for  the  polluted  case  was  offset  by  the  predominantly  forward 


FIG.  19 

FRACTIONAL  DISTRIBUTION  OF  EXTRATERRESTRIAL  INCIDENT 
ENERGY  AT  OMAHA  AS  A  FUNCTION  OF  LOCAL  TIME,  FOR  CLEAR 
AND  INDUSTRIALLY  POLLUTED  ATMOSPHERIC  CONDITIONS 

scattering  by  aerosols.  Evidently,  the  ratio  of  the  diffusely  transmitted 
to  the  directly  transmitted  radiation  woulc  be  increased  by  the  pollutants. 
Absorption  in  the  2,000-foot  polluted  layer  is  a  major  factor  in  the 
polluted  atmosphere,  being  greater  than  the  total  backscatterlng  during 
part  of  the  sunlit  period.  The  actual  absorption  may  have  been  even 
larger,  as  the  assumption  of  independent  and  consecutive  scattering  and 
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absorption  tends  to  minimize  absorption.  The  fraction  transmitted  is 
seen  to  be  always  less  than  85^,  even  for  the  clear  atmosphere. 

Figure  20  is  a  plot  of  the  energy  involved  in  the  various  processes, 
as  a  function  of  local  time.  The  solar  constant  was  taken  as  2.0  ly/min. 
The  amount  of  energy  backscattered  to  space  is  relatively  insensitive  to 
solar  position  except  at  large  solar  zenith  angles,  whereas  there  is  a 
strong  diurnal  variation  of  the  other  components. 


06  06  10  12  14  16  18  20 
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FIG.  20 

DISTRIBUTION  OF  INCIDENT  SOLAR  ENERGY  AT  OMAHA  Ai>  h 
FUNCTION  OF  LOCAL  TIME,  FOR  CLEAR  AND  INDUSTRIALLY 
POLLUTED  ATMOSPHERIC  CONDITIONS 
The  Solar  Constant  was  Taken  as  2.0  ly/min. 


Integration  under  the  curves  of  Fig.  20  yields  the  daily  energy 
distribution  given  in  Table  II  for  the  various  processes.  The  values 
compare  well  with  the  measurements  available  in  similar  conditions. 

From  an  extensive  study  of  daylight  illumination  for  photographic  exposure, 
Jones  and  Condit*®  found  that  light  haze  decreases  the  illuminance  on  a 
horizontal  surface  at  the  ground  by  27^  from  its  value  for  clear  conditions. 
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TABLE  II 


SOLAR  ENERGY  DISTRIBUTION  BY  ATMOSPHERIC  PROCESSES  FOR  THE 
ENTIRE  DAY  FOR  DIFFERENT  ATMOSPHERIC  CONDITIONS  AT  OMAHA 


Atmospheric  Characteristics 

Process 

Energy 

(langleys) 

Percent  of  Incident 
Extraterrestrial 
Flux 

Industrial  pollution  in 

Backscatter 

146.3 

23.8 

lower  2000  feet ,  clear 

Absorption 

131.8 

21.5 

above. 

Transmission 

335.9 

54.7 

Clear 

Backscatter 

133.1 

21.7 

Absorption 

0 

0 

Transmission 

480.9 

78.3 

The  analogous  decrease  here  is  30^.  Another  comparison  can  be  made  by 
combining  the  computed  Rayleigh  backseat  ter,  given  previously  in  Fig.  10, 
with  the  dust  attenuation  factor  of  Houghton,^®  namely,  l-(0.95)"'  where  m 
is  the  air  mass.  Houghton  estimated  that  l/2  of  the  depleted  energy  was 
absorbed,  l/4  backscattered  to  space  and  l/4  diffusely  transmitted. 

These  effects  are  shown  by  the  curves  of  Fig.  21,  Curve  (a)  being-  for 
dust  and  (b)  for  Rayleigh  scattering  for  the  spectral  region 
3100  A  <  X  <  25,000  A.  Curve  (c),  the  sum  of  (a)  and  (b),  represents 
the  total  backscatter  by  the  clear  atmosphere.  This  should  be  compared 
with  Curve  (e),  which  was  transferred  directly  from  Fig.  19.  It  is  seen 
that  only  about  half  of  the  backscatter  computed  from  Waldram's  measure¬ 
ments  is  accounted  for.  In  an  effort  to  obtain  better  agreement,  it  was 
thought  that  Rayleigh  scattering  should  be  represented  by  that  for  only 
the  visible,  since  Waldram's  measurements  were  for  visible  light. 

Curve  (d)  is  the  sum  of  the  fractional  Rayleigh  backscatter  for 
4000  ^  <  X  <  7000  ^  and  the  dust  backscatter  of  Curve  (a).  There  is 
still  considerable  discrepancy  between  Curves  (d)  and  (e),  indicating 
that  the  effects  ofjdust  are  not  adequately  represented  for  this  case  by 
Houghton's  simple  law. 

Blackwell,  et  al.®®*  determined  the  amount  of  absorption  and  scattering 
by  atmospheric  particulates  in  the  clear  atmosphere  at  Kew  Observatory  as 

*  Results  quoted  by  Sheppard®^  are  given  here,  as  the  original  paper  was 
not  available. 
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FIG.  21 

FRACTION  OF  INCIDENT  ENERGY  BACKSCAT7ERED  TO  SPACE 
BY  CLEAR  ATMOSPHERE 


given  in  Table  III.  It  should  be  noted  that  this  is  in  addition  to 
Rayleigh  scattering.  If  these  values  are  added  to  the  Rayleigh  effect, 
given  by  Curve  (c)  of  Fig.  21,  reasonably  good  agreement  with  the  results 
derived  from  Waldram's  measurements.  Curve  (e),  is  obtained. 

TABLE  III 

FRACTIONS  OF  EXTRATERRESTRIAL  VERTICAL  SOLAR  ENERGY  FLUX 
ABSORBED  AND  BACKSCATTERED  BY  ATMOSPHERIC 
PARTICULATES  AT  KEW 


Solar  Zenith  Angle 

36° 

53° 

66° 

•*4 

00 

o 

Fraction  absorbed 

0.06 

0.08 

0.11 

0.16 

Fraction  backscattered  to 

0.05 

0.05 

0.04 

0.02 

space 

Source:  Blackwell,  et.  al.  (Ref.  30). 
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E.  OZONE  ABSORPTION  OF  SOLAR  RADIATION 


A  large  number  of  investigations  of  the  atmospheric  ozone  problem 
have  established  the  main  features  of  the  space  and  time  distribution  of 
ozone  and  its  role  in  the  atmospheric  heat  balance.  The  excellent 
monograph  by  Craig®*  outlines  the  history  of  ozone  work  and  Includes  a 
relatively  complete  bibliography  of  the  field  up  to  1950.  The  main 
developments  since  that  time  have  been  (l)  a  re-determination  of  the 
absorption  coefficients  (Vigroux®®),  (2)  further  measurements  of  the 
amount  and  distribution  of  atmospheric  ozone  as  a  function  of  height, 
geographical  location,  season,  and  weather  situation  (Diitsch,®* 

Ramanathan  and  Kulkarnl,®®  Ramanathan  and  Dave®®),  (3)  more  accurate  data 
on  the  ultraviolet  solar  spectrum  (Johnson®'^),  and  (4)  new  computations 
of  atmospheric  temperature  changes  in  the  ozonosphere  (Murgatroyd  and 
Goody,®®  Pressman®®). 

In  order  to  compute  the  amount  of  energy  absorbed  by  ozone  at  Omaha 
during  the  selected  synoptic  situation,  use  was  made  of  the  3-year  series 
of  umkehr  measurements  at  Arosa  reported  by  Dutsch.®®  The  mean  ozone 
distribution  for  March  was  used  as  being  representative  of  conditions  at 
Omaha  for  the  entire  four-day  period  March  12-15,  up  to  the  limit  of 
Dlitsch's  data  at  47.8  kilometers.  An  exponential  decrease  of  ozone  content 
with  height  was  assumed  from  47.8  km  to  70  km,  the  value  at  70  km  being 
that  given  by  Johnson.®'^  The  ozone  amount  above  70  km  was  assumed  negli¬ 
gible.  Figure  22  shows  the  vertical  distributions  for  Arosa  and 
White  Sands,  New  Mexico,  together  with  the  extrapolation  assumed  for 
Omaha.  The  total  ozone  amount  between  10  and  70  kilometers  was  0.33  cm 
at  NTP. 


The  method  of  Craig®®  was  used  for  computing  the  ozone  absorption 
of  solar  radiation  for  the  sunrise- to-sunset  period  at  Omaha.  Use  of  the 
more  recent  absorption  coefficients  obtained  by  Vigroux®®  would  change 
the  absolute  values  somewhat,  but  not  the  general  patterns  of  ozone 
absorption.  The  total  energy  absorbed  per  unit  volume  at  altitude  z 
is  given  by 


E 


z 


"  j  lox  “x 

J  O 


dX 


(23) 
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FIG.  22 

MEASURED  VERTICAL  DISTRIBUTIONS  OF  OZONE  (SOLID  CURVE) 
AT  AROSA  AND  WHIlt  SANDS,  AND  EXPONENTIAL  EXTRAPOLATION 
(DASHED  CURVE)  ASSUMED  FOR  OMAHA 


where  n  is  the  amount  of  ozone  per  unit  volume,  the  monochromatic 
intensity  at  level  z,  and  the  monochromatic  ozone  absorption  co¬ 
efficient.  The  ozone  pathlength  N  is  given  by  the  relation 


N 


roo 

n  sec  0Q  dz 
J  z 


(24) 


for  the  incident  parallel  solar  radiation.  Effects  of  atmospheric 
scattering  are  neglected  here.  The  rate  of  temperature  change  by  ozone 
absorption  is 


pjj  being  the  air  density  at  level  z,  and  Cp  the  specific  heat  of  air  at 
constant  pressure. 

Curves  of  the  heating  rate  by  ozone  absorption  of  solar  radiation 
as  a  function  of  height  are  shown  in  Fig.  23  for  different  zenith  angles 
of  the  sun.  At  the  highest  levels  the  heating  rate  is  essentially  inde¬ 
pendent  of  solar  position _and  amounts  to  a  few  tenths  of  a  degree  per  hour. 


FIG.  23 

RATI  OF  HEATING  BY  OZONE  ABSORPTION 
OF  SOLAR  RADIATION  AS  A  FUNCTION  OF 
HEIGHT  FOR  DIFFERENT  SOLAR  ZENITH  ANGLES 

The  maximum  rate  is  seen  to  be  1.7*C/hr  for  zenith  sun.  The  noon  solar 
zenith  angle  in  March  at  Omaha  is  about  47“ ,  for  which  the  maximum  heating 
rate  of  1.5“c/hr  occurs  at  47  km. 

A  time  integration  of  the  heating  rate  for  the  sunrise-to-sunset 
period  at  Omaha  results  in  the  diurnal  temperature  increase  by  ozone 
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absorption  as  a  function  of  height  given  in  Fig.  .'24.  The  12.5*C  maximum 
is  considerably  greater  than  the  8®C  maximum  given  by  Murgatroyd  and 
Goody,®*  reflecting  the  greater  total  ozone  content  assumed  here.  The 
50  km  height  of  the  maximum  agrees  well  with  their  results,  and  with  the 
well  known  temperature  maximum  in  that  region. 


HEATING  INCREMENT - ‘C 

FIG.  24 

DIURNAL  TEMPERATURE  INCREASE  BY  OZONE  ABSORPTION 
OF  SOLAR  RADIATION  AT  OMAHA  IN  MARCH  1957, 

AS  A  FUNCTION  OF  HEIGHT 

The  rate  of  energy  absorption  by  the  total  ozone  content  of  the 
atmosphere  at  Omaha  for  March  15  is  shown  as  a  function  of  local  time  by 
the  curve  of  Fig.  25.  The  maximum  of  0.285  ly/min  occurring  at  noon 
represents  2.1^  of  the  extraterrestrial  incident  energy.  Integration 
over  the  sunlit  period  shows  a  total  absorption  by  ozone  of  13.55  ly, 
which  is  2.2^  of  the  total  incident  solar  energy  for  the  period. 
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FIG.  25 

RATE  OF  SOLAR  ENERGY  ABSORPTION  BY  OZONE  AT  OMAHA 
IN  MARCH  1957,  AS  A  FUNCTION  OF  LOCAL  TIME 
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APPENDIX 
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THE  MONOCHROMATIC  REFLECTANCE  Ry  OF  THE  ATMOSPHERE-SURFACE  SYSTEM 
FOR  VARIOUS  VALUES  OF  THE  OPTICAL  THICKNESS  t  AND 
SURFACE  REFLECTIVITY,  A 


m 

x(X) 

Surface  Reflectivity,  A 

D 

1^0 

0 

0.05 

0.10 

0.15 

0.25 

0.50 

0.80 

1.0 

3120 

0.02 

0.751 

0.758 

0.773 

0.790 

0.840 

0.922 

0.10 

0.704 

0.712 

HI 

0.730 

0.750 

0.809 

0.908 

0.650 

0.660 

0.681 

0.705 

0.775 

0.891 

0.40 

0.546 

0.558 

0.572 

0.586 

0.616 

0.707 

0.859 

0.60 

0.458 

0.474 

0.490 

0.507 

0.543 

0.651 

0.831 

0.80 

0.409 

0.427 

0.446 

0.486 

0.608 

0.811 

0.92 

0.359 

0.377 

0.396 

0.416 

0.469 

0.587 

0.801 

1.00 

0.340 

0.359 

0.378 

0.399 

0.443 

0.575 

0.795 

0.50 

3660 

0.02 

0.675 

0.687 

0.698 

0.915 

0.10 

0.610 

0.624 

0.638 

0.653 

Ha 

0.898 

0.20 

0.523 

0.540 

0.558 

0.576 

0.614 

Ela 

0.875 

0.40 

0.383 

0.405 

0.428 

0.451 

0.500 

0.638 

0.838 

0.60 

0.297 

0.322 

0.348 

0.374 

0.430 

0.587 

0.816 

0.80 

0.241 

0.268 

0.296 

0.325 

0.385 

0.554 

0.801 

0.92 

0.216 

0.244 

0.273 

0.303 

0.365 

0.540 

0.795 

1.00 

0.202 

0.231 

0.260 

0.291 

0.354 

0.532 

0,791 

0.25 

4365 

0.02 

0.609 

0.625 

0.641 

0.658 

0.785 

0.909 

0.10 

0.503 

0.524 

0.545 

0.566 

0.610 

0.727 

0.804 

0.20 

0.375 

0.401 

0.427 

0.454 

0.509 

0.657 

0.854 

0.40 

0.238 

0.270 

0.302 

0.335 

0.402 

0.582 

0.822 

0.60 

0.173 

0.207 

0.242 

0.278 

0.351 

0.546 

0.807 

0.80 

0.136 

0.172 

0.208 

0.245 

0.321 

0.525 

0.798 

0.92 

0.120 

0.157 

0.194 

0.231 

0.309 

0.517 

0.794 

1.00 

0.112 

0.148 

0.186 

0.224 

0.302 

0.512 

0.792 

0.15 

4950 

0.588 

0.667 

mmi 

0.406 

0.432 

0.459 

0.486 

0.541 

0.869 

iIETim 

0.270 

0.302 

0.335 

0.368 

0.436 

0.612 

0.839 

0.158 

0.195 

0.233 

0.271 

0.349 

0.552 

0.814 

0.60 

0.111 

0.151 

0.191 

0.231 

0.313 

0.528 

0.803 

0.80 

0.086 

0.126 

0.167 

0.209 

0.293 

0.514 

0.798 

0.92 

0.076 

0.116 

0.158 

0.200 

0.285 

0.508 

0.796 

1.00 

0.070 

0.111 

0.153 

0.195 

0.281 

0.505 

0 . 794 

Note;  The  solar  zenith  angle  0q  is  given  in  terms  [Iq  =  cos  0q. 
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6440  0.02 


.02 

.10 

.20 

.40 

0.60 

0.80 

0.92 

1.00 


0.012 

0.008 

0.006 

0.005 

0.005 


0.006 

0.002 

0.001 

0.001 

0.000 


Surface  Reflectivity,  A 


0.05  0.10  0.15  I  0.25  0.50  I  0.80 


.541 
.323 
.199 
0.111 
0.60  0.077 

0.80  0.059 

0.92  0.052 

1.00  0.048 


0.562 

0.355 

0.236 

0.152 

0.119 

0.102 

0.095 

0.091 


.583  0.605 

.386  0.418 

.273  0.310 

.193  0.235 

.162  0.205 

,146  0.190 

139  0.183 

136  0.180 


_  0.760  0 

.482  0.647  0 

.386  0.582  0 

.319  0.536  0 

.293  0.518 

.279  0.509 

.273  0.505 

.270  0.503 


.902 

.855 

.828 

0.809 

0.802 

0.798 

0.797 

0.796 


0.419  0.486  0.656 


C 

0.225  0.315  I  0.541 

0.094  I  0.141  I  0.188  0.282  0.519 


.40  0.059 

0.60  0.040 

0.80  0.030 

0.92  0.026 

1.00  0.024 


0.072 

0.065 

0.061 

0.059 

0.059 


.352 
.136 
0.094 
0.072 
0.065 
0.061 
0.059 
0.059 


.237 
.095 
.073 
0.061 
0.057 
0.055 
0.055 
0.054 


0.115 
0.Q63 
0.056 
0.052  I  0.102 


0.120  0.168 
0.113  0.162 

0.110  0.158 

0.108  0.157 

0.107  0.156 


0.265  0.508 
0.259  0.504 
0.256  0.501 
0.255  0.501 
0.254  0.500 


0.141 


0.120  0.168 
0.113  0.162 

0.110  0.158 

0.108  0.157 

0.107  0.156 


.419  0.486 

.225  0.315 

.188  0.282 
0.265 
0.259 
0.256 
0.255 
0.254 


0.656 

0.541 

0.519 

0.508 

0.504 

0.501 

0.501 


0.862 
0.815 
0 . 807 
0.802 
0.800 
0.799 
0.799 


0.500  0.799 


II 


.110  0.159 

0.107  0.156 


.170  0.267 

.159  0.257 


0.254 

0.253 


0.051 

0.050 

0.050 


0.105  0.154  0.253 

0.104  0.153  0.252 

0.104  0.153  0.252 


0.208 

0.161 

0.155 

0.152 

0.151 


0.597 

0.521 

0.510 

0.504 

0.502 

0.501 

0.500 


0.838 

0.808 

0.803 

0.801 

0.800 

0.800 

0.800 


0.500  0.799 


0.534 


0.506 

0.502 


2 

0.801 


0.101 


0.251 

0.250 


0.500  0.800 

0.500  0.800 


0.100  0.150  0.250  0.500 

0.100  0.150  0.250  0.500 


0.000  0.050  0.100  0.150  0.250 


0.800 
0.800 
0.500  I  0.800 


Appendix  (Concluded) 


/ 


T 

x(A) 

Ho 

Surface  Reflectivity, 

A 

0.05 

0.10 

0.15 

0.25 

0.50 

0.80 

0.001 

17000 

Esa 

0.023 

0.121 

.0.512 

0.805 

EBl 

0.004 

0.054 

0.104 

0.502 

0.801 

0.002 

0.101 

rgm 

0.251 

0.501 

0.800 

0.40 

0.000 

0.050 

0.100 

0.150 

0.250 

0.500 

0.800 

0.60 

0.000 

0.050 

0.100 

0.150 

0.250 

0.500 

0.800 

0.80 

0.000 

0.050 

0.100 

0.150 

0.250 

0.500 

0.800 

0.92 

0.000 

0.050 

0.100 

0.150 

0.250 

0.500 

0.800 

1.00 

0.000 

0.050 

0.100 

0.150 

0.250 

0.500 

0.800 
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